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Abstract. We study limit joint distributions of the ensemble of symmetric blocks of 
independent Gaussian random matrices with block-identical variances (Gaussian Sym- 
metric Block Ensemble). Our approach is based on the concept of noncommutative 
independence called matricial freeness which reminds freeness in free probability, but 
describes the asymptotics of symmetric blocks rather than that of the whole matri- 
ces. The main result can be viewed as a block refinement of Voiculescu's asymptotic 
freeness of the ensemble of independent Gaussian random matrices. We show that 
this concept provides a unified framework for studying the asymptotic distributions 
of sums and products of independent Gaussian random matrices, including random 
matrices of Wishart type and random matrix models for free Bessel laws. Moreover, 
it leads to random matrix models for boolean independence, monotone independence 
and s-freeness. 



1. Introduction and main results 

One of the most important features of free probability is its close relation to ran- 
dom matrices. It has been shown by Voiculescu that independent Hermitian Gaussian 
random matrices are asymptotically free [32]. This result shows that the concept of 
freeness, or free independence, is not only a noncommutative analog of classical inde- 
pendence, but is also fundamental in the random matrix theory, which puts the classical 
result of Wigner on the semicircle law as the limit distribution of certain symmetric 
random matrices [36] in an entirely new perspective. Moreover, it is universal in view 
of Dykema's generalization to non-Gaussian random matrices [12]. 

In particular, if we are given an ensemble of independent Hermitian n x n random 
matrices 

{Y(u, n) : ue U) 

whose entries are suitably normalized and independent complex Gaussian random vari- 
ables for each natural n, then 

lim r(n)(F(wi,n) . . .Y(u m ,n)) = $(u(ux) ■ ■ -co(u m )) 

n— >ao 

for any u\, . . . , u m e U, where {co(u) : u e U) is a semicircular family of free Gaussian 
operators living in the free Fock space with the vacuum state <£> and r(n) is the normal- 
ized trace composed with classical expectation. This realization of the limit distribution 
gives a fundamental relation between random matrices and operator algebras. 
Consequently, if Zl is finite, the moments of the sum 

Y(n) = YjY(u,n) 
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under r(n) converge to the moments of the sum of free Gaussian operators 

■iieu 

and can be described in terms of the free additive convolution of semicircle laws 

a = o-(l)fflo-(2)EB...EBo-(m) 

where m is the cardinality of 11. The free additive convolution, introduced by Voiculescu 
[33] for compactly supported measures on the real line and generalized by Bercovici 
and Voiculescu [7] to measures with non-compact support, turned out to be an effective 
tool in the study of asymptotic spectral distributions of sums of independent Hermitian 
random matrices. 

The basic original random matrix model studied by Voiculescu corresponds to inde- 
pendent complex Gaussian variables, where the entries Y it j(u, n) of each matrix Y(u, n) 
satisfy the Hermiticity condition Y it j(u,n) = Yj ti (u,n) as well as 

E(Yij(u,n)) = and EQY^u, n)\ 2 ) = l/n 

for any 1 ^ i,j ^ n and u e 11. If we relax the second assumption, the scalar- valued 
free probability is no longer sufficient to describe the asymptotics of Gaussian random 
matrices. One approach is to use the operator- valued free probability, as in the work of 
Shlakhtyenko [30] , who studied the asymptotics of Gaussian random band matrices and 
proved that they are asymptotically free with amalgamation over some commutative 
algebra. 

Our approach is based on the decomposition of independent Hermitian Gaussian 
random matrices Y(u,n) with block-identical variances of \Y it j(u,n)\ into symmetric 
blocks T p q (u, n), where ue 11 and n e N, namely 

Y(u,n) = Yj T P,q( U i n ) 

where each T PA {u,n) is built from the entries Y it j(u,n) indexed by pairs taken 
from the union (N p x iV 9 ) u (N q x N p ), where 

\N | 

[n] = Ni u . . . u N r , with lim — — = d q ^ for any q, 

n^oo n 

is the partition of the set [n] := {1,2, . . . ,n} into disjoint non-empty subsets (the fact 
that they depend on n is supressed in the notation). 

At the same time, we decompose each semicircular Gaussian operator in terms of the 
corresponding r x r array of matricially free Gaussian operators [23], namely 

r 

w(w) = J] Wp, g («), 
p,q=l 

living in the matricially free Fock space of tracial type, similar to the matricially free 
Fock space introduced in [23], in which we distinguish a family of states . . . ,^ r }, 
used to build the array (& p , q ) by setting ty Ptq = ^/ q . 

Nevertheless, in order to reproduce the limit distributions of symmetric blocks in the 
case when the variances of \Y it j(u, n)\ are only block-identical rather than identical, we 
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still need to rescale the matricially free Gaussian operators in the above expressions. 
The corresponding arrays of distributions 



[a(u)] 



( 0-1,1 («) K\,l{u) ■ ■ ■ Kl,r{u) \ 
K2,l(u) CT 2i 2(m) • • • K 2 , r (u) 

\ Kr,i(u) Kt,2(m) • • • Vr,r{u) J 



consisting of semicircle laws a q>q (u) and Bernoulli laws K, p>q (u), replace the semicircle 
laws <t(u) and play the role of matricial semicircle laws. In turn, the rescaled matricially 
free Gaussians (u] Piq (u)) replace free Gaussians u(u) and thus the corresponding sums 
become more general objects called random pseudomatrices. 

In any case, it is the ensemble of symmetrized Gaussian operators defined by 



Uq,q(u) if p = q 

co p , q (u) + u} q:P (u) if p ^ q 



which gives the operatorial realizations of the limit joint distributions of the ensemble 
of symmetric blocks of the family {Y(u,n) : u s 11} . The array of distributions cor- 
responding to each (Q Piq (u)) consists of diagonal semicircle laws and off-diagonal laws 
defined by two-periodic Jacobi sequences which becomes semicircle laws if the second 
moments of cu P;q (u) and u q>p {u) are equal. 

In this general setting, we prove that the ensemble of symmetric blocks of independent 
Hermitian Gaussian random matrices 

{T, m (u,n) :p,qe[r],ueU} 

converges in moments as n — > go to the ensemble of symmetrized Gaussian operators 
under the normalized partial traces, namely 

lim T q (n)(T Pim (u 1 ,n)...T Pmtqm (u m ,n)) = ^ q {ui Puqi { Ul ) . . .u PmSm {u m )), 

n— >oo 

where u\, . . . , u m e Zl, all remaining indices belong to [r] := {1,2,..., r}, and T q (n) 
denotes the normalized partial trace over the set of basis vectors {e^ : k e iVq} composed 
with classical expectation (called a partial trace in the sequel). 

This result can be viewed as a block refinement of that used by Voiculescu in his 
fundamental asymptotic freeness result [32]. Using the convex linear combination 

r 

* = 2 d q %, 

5=1 

we easily obtain a similar formula for the normalized trace composed with classical 
expectation r{n) (called the trace in the sequel) considered by most authors in their 
studies of asymptotics of random matrices. It is obvious that asymptotic freeness is a 
special case of our asymptotic matricial freeness and it corresponds to the case when 
the variances of all |Yy(-u,n)| are identical. However, we also show that by consid- 
ering partial traces we can produce random matrix models for boolean independence, 
monotone independence and s-freeness. It is not a coincidence since all these notions of 
independence arise in the context of suitable decompositions of free random variables 
as shown in [21,22]. 
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Using the concept of matricial freeness, we studied asymptotic distributions of sym- 
metric blocks of one Hermitian Gaussian random matrix in [23]. In this paper, we gen- 
eralize that result in two important directions. Firstly we consider symmetric blocks of 
the ensemble of independent Hermitian and non- Hermitian Gaussian random matrices. 
Secondly we include the case when some elements of the dimension matrix 

D = diag(di, d 2 , ■ ■ ■ , d r ), 

called asymptotic dimensions, are equal to zero. This corresponds to the situation when 
at least one dimension of some blocks does not grow fast enough as n — > go. Of special 
interest are off-diagonal symmetric blocks consisting of two rectangular blocks whose 
dimensions do not grow proportionately with the larger dimension proportional to n and 
the smaller dimension growing slower than n (such blocks will be called unbalanced). 
In the formulas for joint distributions under the trace r(n), the contributions from 
partial traces associated with vanishing asymptotic dimensions disappear in the limit. 
However, certain limit joint distributions involving unbalanced symmetric blocks under 
the partial traces themselves become non-trivial and interesting. 

Limit distributions of the sums Y(n) of independent Hermitian Gaussian random 
matrices can be described in terms of convolutions of arrays of matricial semicircle 
laws. Namely, one can define the convolution of arrays of distributions, where the 
diagonal distributions and the off-diagonal ones are convolved by means of free and 
boolean convolutions, respectively. Thus 

[a] = [^(l)] ffl ^(2)] B ... H [(7(m)] 

is a convolution of matricial semicircle laws, where the notation EH is used since this 
convolution is a generalization of the free additive convolution to the category of arrays 
of distributions. In turn, the limit distributions of Y(n) under the trace are obtained 
by taking the matricially free convolution [*] of the entries of [a] . This convolution is 
introduced and studied in this paper. 

The convolution EB is related to the strongly matricially free convolution introduced 
and studied in [24] . We can say that it describes the asymptotic relation between sym- 
metric blocks having the same position in the considered independent matrices Y(u, n), 
whereas the convolution [*] describes the asymptotic relation between symmetric blocks 
having different positions in the matrix Y{n). Putting two convolutions together, we 
obtain a formula called the Double Convolution Formula. 

Next, we derive the limit joint distributions of the ensemble of symmetric blocks of 
independent non-Hermitian Gaussian random matrices as it was done by Voiculescu 
in the case of independent non-Hermitian Gaussian random matrices themselves [32]. 
Further, by considering off-diagonal symmetric blocks, we can recover limit joint distri- 
butions of independent Wishart matrices 

W(u,n) = B(u,n)B*(u,n) 

where each B(u,n) is a Gaussian random matrix [37], as well as distributions of matrices 
of similar form W(n) = B(n)B*(n), where B{n) is a sum or a product of independent 
Gaussian random matrices. In all these situations, it suffices to embed the considered 
matrices in the algebra of symmetric blocks in an appropriate way. In particular, we 
obtain the random matrix model for the noncommutative Bessel laws of Banica et al [4]. 
Non-Hermitian Wishart matrices, studied recently by Kanzieper and Singh [19], can be 
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treated in a similar way. We also expect that a continuous generalization of our block 
method should allow us to treat symmetric triangular random matrices studied by Basu 
et al [5] and the triangular ones studied by Dykema and Haagerup [13]. An extension of 
our results to the class of non-Gaussian random matrices (under suitable assumptions, 
like uniform boundedness of type given in [12]) seems to present no difficulty 

The paper is organized as follows. Section 2 is devoted to the concept of matricial 
freeness, the corresponding arrays of matricially free Gaussian operators and matricial 
semicircle laws. Symmetrized counterparts of these objects are discussed in Section 3. 
In Section 4, we describe the combinatorics of the mixed moments of operators of both 
types. The first main result of this paper is contained in Section 5, where we derive the 
limit joint distributions of the Gaussian Symmetric Block Ensemble under partial traces. 
The corresponding convolution limit theorem with the Double Convolution Formula is 
proved in Section 6. In Section 7, we apply this formalism to construct random matrix 
models for asymptotic boolean independence, monotone independence and s-freeness. 
In Section 8, we study the transforms of limit distributions and introduce matrix- valued 
analogs of Catalan numbers. Non-Hermitian Gaussian random matrices are treated in 
Section 9. Section 10 is devoted to Wishart matrices, matrices of Wishart type and 
products of independent Gaussian random matrices. 

Finally, we take this opportunity to issue an erratum to the definition of the symmet- 
rically matricially free array of units given in [23]. The new definition entails certain 
changes which are discussed in the Appendix. 

2. Matricial semicircle laws 

Let us recall the basic notions related to the concept of matricial freeness [22]. Let A 
be a unital *-algebra with an array (Aij) of non-unital *-subalgebras of A and let (<fi,j) 
be an array of states on A. Here, and in the sequel, we shall skip J in the notations 
involving arrays and we shall tacitly assume that e J, where 

j c [ r ] x [r], 

and r is a natural number. Further, we assume that each Aij has an internal unit 
by which we understand a projection for which alij = l^a = a for any a e Aij, and 
that the unital subalgebra X of A generated by all internal units is commutative. 

For each natural number m, let us distinguish the subset of the m-fold Cartesian 
product J x J x ... x J of the form 

5m = • • • , (im,jm)) ■ # • • • # («m, jm)}, 

and its subset 

%.m = {((iuh), (i2,*3) ■ ■ ■ , (imjm)) ■ # fah) # ■ ■ ■ # (*m, jm)}- 

In other words, the neighboring pairs of indices in the set 1<i m are not only different (as 
in free products), but are related to each other as in matrix multiplication. Objects 
labelled by diagonal or off-diagonal pairs, respectively, will be called diagonal and off- 
diagonal. 

If fi, y?2, ■ ■ ■ ,(p r are states on A, we form an array of states ((pij) as follows: 
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and then we will say that ((fiij) is defined by the family (<Pj). In particular, when 
3 = [r] x [r], this array takes the form 



( Ifii tp 2 ... ip r \ 
<Pl </?2 ••• <Pr 



\ (pi ip 2 ... (fir J 

The states tpi, ip 2 , ■ ■ ■ <p r can be defined as states which are conjugate to a distinguished 
state (p on A, where conjugation is implemented by certain elements of the diagonal 
subalgebras, namely 

(Pj(a) = fitfabj) 

for any a e A, where bj e Ajj n Keiip is such that ip(b*bj) = 1 for any j e [r]. 

In general, other arrays can also be used in the definition of matricial freeness, but in 
this paper we will only use those defined by a family of r states, say (<fj), and thus the 
definition of matricial freeness is adapted to this situation. Since the original definition 
uses the diagonal states [22], which here coincide with (fi, . . . , <p r , we shall use the latter 
in all definitions. 

Definition 2.1. We say that (ljj) is a matricially free array of units associated with 
(Aij) and (tpij) if for any state tpj it holds that 

(a) ipj(biab 2 ) = <Pj(bi)(pj(a)ipj(b 2 ) for any as A and bi, b 2 e X, 

(b) <Pj(l k ,i) = for any i,j, k, I, 

(c) if a r e A ir j r n Kerip irjjr , where 1 < r < m, then 

,„ („ 1 „ „ \ _ f ^j(« a 2 • • -dm) if ((«1, Jl), • • • , (imJm)) e ^C m 

^al n) ,a 2 ...a mj -| otherwise 
where a e ^4 is arbitrary and ((ii, Ji), - - - , (i m ,jm)) e J7 m - 

Definition 2.2. We say that *-subalgebras (^4jj) are matricially free with respect to 
(<ftj) if the array of internal units (ljj) is the associated matricially free array of units 
and 

ipj{a x a 2 . . . a n ) = whenever a k e A ikJk n KerLp ikJk 

for any state ipj, where • • • , (i m ,jm)) s 3m- The matricially free array of vari- 

ables (ctij) in a unital *-algebra A is defined in a natural way. 

The most important example of matrically free operators are the matricial gener- 
alizations of free Gaussian operators. They are defined on a suitable Hilbert space of 
Fock type. In order to define this Fock space, we shall use boolean and free Fock spaces. 
Recall that by the boolean and free Fock spaces over the Hilbert space H, respectively, 
we understand the direct sums 

00 

T (H) = CVl®n and T(U) = rae© u® m , 

ra—1 

where Q is a unit vector. Both spaces are endowed with the canonical inner products. 

Let us assume now that J7 contains the diagonal and that 11 is a finite index set 
(in many instances, these assumptions can be dropped). Moreover, in the operatorial 
context discussed below, we shall use indices (p, q) instead of since the related 
operators give the limit realizations of symmetric blocks T p q (u, n) built from complex 
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Gaussian variables Yij(u,n) indexed by pairs falling into the union of two Carte- 
sian products of intervals N p and N q . 

Thus, to each (p, q) e J and u e U we now associate a Fock space according to the 
following rule: 

T(Ce g>g (u)) if p = q 
To(Ce P!q (u)) if p ^ q ' 

where e Ptq (u) is a unit vector for any (p, q) e J and ue U. Each Fock space Fp^u) has 
a distinguished unit (vacuum) vector denoted Q Ptq (u). In other words, we put free Fock 
spaces over one- dimensional Hilbert spaces on the diagonal and boolean Fock spaces 
(over one-dimensional Hilbert spaces) elsewhere. This special type of assignment, which 
appeared in our previous work [23], is related to the difference between multiplication 
of diagonal and off-diagonal entries of random matrices. 

Definition 2.3. By the matricially free Fock space of tracial type we will understand 
the direct sum of Hilbert spaces 

M = ®M q , 

9=1 

where each summand is of the form 



• (Ur 



Mq = cn q ®(& © ® • • • ® K 

where J- P:q {u) is the orthocomplement of CQ Ptq (u) in Tp^u), endowed with the canon- 
ical inner products. 

Each space A4 q is similar to the matricially free-boolean Fock space introduced in 
[23] except that its tensor products do not have to end with diagonal free Fock spaces 
and that we use a family of arrays of Fock spaces 

as well as a family of vacuum vectors 

instead of one array and one distinguished vacuum vector. Of course, each Q Ptq (u) is 
different from any vacuum vector from the set Q. The states defined by . . . , Q r will 
be denoted ^/ r , respectively. Thus, 

for any a 6 B{M) and any q e [r]. By (& Pt q) we shall denote the array defined by the 
family thus = ^ q for any (p, q) e J. 

There exists a canonical embedding of M in the free Fock space over the Hilbert 
space with the orthonormal basis {e p>g (u) : (p,q) e J,u e 1/}, except that the single 
vacuum vector in this free Fock space is replaced by the family Q, namely 

r:M- 0CO, Ce p , q (u)Y 

\q=l J \(p,q,u)ejxzi J 

These embeddings will be used in the definition of our matricially free Gaussian oper- 
ators which remind free Gaussian operators, but they have larger kernels. 
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Remark 2.1. The orthonormal basis B of the matricially free Fock space M is given 
by vacuum vectors fl\, . . . , Q r , and simple tensors of the form 

<&%M ® eg5,(«2) ® ■ ■ ■ ® 
where (pi,P2), ■ • • , (p m) e .7, Ui, ■ ■ ■ , u m e U and l lr ..,i m ,meN, for which it holds 
that 

(Pl,P2,«l) ^ (P2,P3, "2) # • • • # (Pm,5,«m), 

with the additional condition that e N or Ik = 1 if is diagonal or off- 

diagonal, respectively, for any k. In particular, the above conditions imply that we 
allow the triples (p,p,u) and (p,p,u f ) to be neighbors if m ^ u', as well as triples 
(p,j,u) and (j,k,u) to be neighbors as long as (p,j) ¥^ (j,k). 

Definition 2.4. Let (otij(u)) be an array of non-negative real numbers for any u e 11. 

We associate with each such matrix the matricially free creation operators 

where (p,q) e J7 and u e 11, and where £(e P:q (u))Q s = 8 q , s e Ptq (u), otherwise £(e Pjq (u)) 
agrees with the free creation operator associated with the vector e Ptg (u). The corre- 
sponding sums 

= PpM + pp, q ( u ), 

where p Ptq (u) is the adjoint of p p , q (u) (called a matricially free annihilation operator), 
will be called matricially free Gaussian operators. The sum 

u(u) = ^^pM 
p,q 

will be called a Gaussian pseudomatrix, with [u;(u)] denoting the matrix (co Ptq (u)). 

Note that the above definition is very similar to that in [23] and for that reason we 
use similar notations and the same terminology. Strictly speaking, however, the above 
operators live in a larger Fock space, built by the action of a collection of arrays of 
creation operators, namely 

{PpM (P,Q) e J,ueU} 
onto the family of vacuum vectors Q rather than by the action of one array onto one 
distinguished vacuum vector. Besides, in contrast to [23], we allow a p>q (u) to vanish, 
which leads to trivial operators. 

Proposition 2.1. The action of p Ptq (u) onto the basis B is given by 

P P , q (u)tt q = a p , q (u)e p , q (u) 

P P ,?(w)(e g , a (*) ® w ) = ct P ,q( u )( e P,q( u ) ® e,, a (t) ® w) 
for any e qjS (t) ®w<eB, where (p, q), (q, s) e J and u,t e U, and in the remaining cases 
the action gives zero. 

Proof. These formulas are easy consequences of Definition 2.4. ■ 

In other words, up to a scalar multiplication, p p , q {u) is given by tensoring on the 
left by e Pjq (u) restricted to the subspace of M. which is matricially related to (p, q). 
The commutation relations for the matricially free creation and annihilation operators 
generalize those for their free counterparts as we show below. 
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Example 2.1. If J consists of one element and a ltl (u) = 1 for any u e U, we can write 
Pi,i( u ) = ^(u), P*i( u ) = and fix = Q. Then we obtain the action of free creation 

operators onto the canonical basis of the free Fock space 

£(u)Q = e(u) and £(u)(e(t) ® w) = e(u) e(t) ® w 

and the well-known relation between free creation and annihilation operators 

t(u)£(t) = 5 u>t 

for any u,t e U, of which the matricially free creation and annihilation operators are 
natural generalizations. 

Example 2.2. The moments of each operator u p>q (u) will be easily obtained from 
Proposition 2.3 (cf. [23, Proposition 4.2]). Therefore, let us compute some simple 
mixed moments. First, we take only off-diagonal operators: 

= t>l p (t)b Ptq (u) 

where p q k ^ p and t,u e U are arbitrary and where b Pj0/ (u) = a pq (u). If some 
operators are diagonal, we usually obtain more terms. For instance 

^ q(^ q (u)u^ p (t)uj ptq (u)) = * 9 (p*,(u)p* i p(t)pp 1 p(t)p* i p(t)pp 1 p(t)pp,g(M)) 

+ ^g(pp, g ( w )^P,p( t )^,p( t )fep( t )fep( t )fe'?( w )) 
= Kp^pM 

for any p ^ q and any u, t e U. Let us observe that the pairings between the diagonal 
annihilation and creation operators and between the off-diagonal ones are the same as 
those between annihilation and creation operators in the free case and in the boolean 
case, respectively. 

Now, using the fact that U is finite, we would like to prove matricial freeness of 
collective Gaussian operators of the form 

Up,q = Pp,q + P p , q , 

where 

ueu ueu 

with respect to the array (^ p>q )- As in [23, Proposition 4.2], we will prove a more 
general result on matricial freeness of the array of algebras (A Ptq ) generated by collective 
creation, annihilation and unit operators. 

For that purpose, let us define a suitable array of collective units (l p , g ). The easiest 
way to do it is to say that is the orthogonal projection onto the subspace onto which 
the associated collective creation or annihilation operators act non-trivially. However, 
a more explicit definition will be helpful. 

Definition 2.5. By the array of collective units we will understand (l p , g ), where 

lp, q = r p , q + s p>q for (p,q) 6 J, 
where s Pjq is the orthogonal projection onto 

Np, q ■= ®Ran(p M («)), 

ueu 
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and r PA is the orthogonal projection onto 

K q , q = Cfi,©@^, 

k^q 
k 

for any diagonal (q, q) e J and off-diagonal (p, q) e J, respectively. 
Proposition 2.2. The following relations hold: 

p p , q (u)p k! i(t) = and p^ q (u)p* kl (t) = 
for q k and any u,t e 11. Moreover, 

■4>w«) - { £fe %l\ 

and otherwise p* q (u)pk,i(t) = for any (p,q,u) ¥^ (k,l,t) e J x 11, where b Piq (u) = 

Proof. These relations follow from Proposition 2.1. ■ 

By the semicircle law of radius 2a > we understand the continuous distribution on 
the interval [—2a, 2a] with density 

da = V '4a 2 — x 2 /2irdx 

and the Cauchy transform 

z- \Jz 2 - Aa 2 

G{Z) = 2a~ 2 ' 

where the branch of the square root is chosen so that \l z 2 — Aa 2 > if z e R and 
z e (2a, go). In turn, the discrete distribution 

k = 1/2 (<L Q + £ a ) 

with the Cauchy transform 

G(z) ' 



z — a 2 1 z 

defines the Bernoulli law concentrated at +a. 
Proposition 2.3. For finite 11 and any (p,q) e J, let 

and let (^ p , q ) be the array of states on B(M) defined by the states ^i, . . . , \& r associated 
with the vacuum vectors Qi, . . . ,Q r , respectively. Then 

(1) the array (A p>q ) is matricially free with respect to (^ p , q ), where the array of units 
is given by Definition 2.5, 

(2) the ^ q>q - distribution of non-trivial ui q>q {u) is the semicircle law of radius 2a q>q {u) 
for any (q, q) e J and us 11, 

(3) the ^ pq - distribution of non-trivial uj pa (u) is the Bernoulli law concentrated at 
+a Ptg (u) for any off-diagonal (p, q) e J and ue 11. 
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Proof. We shall prove only (1) since (2) and (3) follow from the case of one array in 
[23, Proposition 4.2]. For an off-diagonal pair (p,q) e J, we have 

uell 

where b PiQ = Yj u b p , q ( u )- 111 turn, for a diagonal pair (q, q) e J, we have 

ueti 

where b qs = Yjueu^i,<i( u )- The ^ rs ^ relation implies that if p P)Q is non-trivial, then 
r p,qi s p,q e ^p,q whenever p q- Thus, any element a e A pq := ,4p j9 n Ker^^ is of the 
form 

« = ap P , g + f3pp, q + 7 s p, 9 

for some constants a,j3,"f. In turn, if <p q ^ q is non-trivial, any element a e A qq : = 
.4.^^ n Ker^ ?)? is spanned by l g ,q and products of the form 

p?>0 ■ ■ ■ pZ( u M7( v ^) ■ ■ ■ pJ? («;+*), 

where ^ n « + 2; m ? > an d Ui # • • • ^ ^j+fe (one uses commutation relations to obtain 
this form as in [23, Proposition 4.2]). Both facts allow us to conclude that any element 
of A q q , when acting onto a simple tensor w, either gives zero or a linear combination 
of vectors of the form 

if w begins with some e g , t (it), whereas any element of the off-diagonal A pq , when acting 
onto a simple tensor w which begins with some e qtt (u), gives a linear combination of 
vectors of the form e Piq (t)®w for some t and otherwise gives zero. Therefore, a succesive 
action of elements cii, . . . ,a n satisfying the assumptions of Definition 2.2 either gives 
zero or we have 

ai . . . a n Q q £ T p \ m ® . . . ® T v \ q 1 CQ q , 

which implies that 

* 9 (ai ...a n ) = 

for any q e J, which gives the freeness condition of Definition 2.2. The unit conditions 
of Definition 2.1 are checked by direct computation. In particular, \ pq leaves invariant 
any vector from J® 1 ® . . . (g) whenever q = pi and kills it otherwise. This proves 
condition (3) of Definition 2.1. Checking the remaining conditions is straightforward. 
This completes the proof. ■ 



Corollary 2.1. If J = [r] x [r], then the array of distributions of non-trivial u PA (u) , 
us U, in the states {^ p , q ) takes the matrix form 

/ <7i,i (tt) Ki t2 (u) . . . K 1>r (u) \ 

k 2 ,i(u) a 2 , 2 {u) ... K2, r (u) 

K«)] = 

V «r,2(«) • • • <*rA U ) / 

where cr m (u) is the semicircle law of radius 2a qA {u) for any q and K p , q ( u ) ^ s the 
Bernoulli law concentrated at +a p>q {u) for any p ^ q. 
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Proof. This is a special case of Proposition 2.3. ■ 

Each matrix of the above form plays the role of a matricial semicircle law which is 
standard if each diagonal law is the semicircle law of radius two and each off-diagonal 
law is the Bernoulli law concentrated at +1. In fact, the family {[cr(u)],u e U] can 
be treated as a family of independent matricial semicircle laws. This aspect will be 
discussed in Section 6. In the case when some of the operators u) p>q (u) are equal to zero, 
the corresponding measures are replaced by 5o. 

3. Symmetrized matricial semicircle laws 

We will show in Section 5 that the limit distributions of symmetric blocks of Hermit- 
ian Gaussian random matrices can be expressed in terms of matricially free Gaussian 
operators in the most general case studied in this paper, including the situation when 
the dimension matrix D is singular. Nevertheless, if both asymptotic dimensions asso- 
ciated with a given symmetric block T Ptq (u, n) are non-zero, its limit realization is not 
given by u) Ptq (u), but by its symmetrization uj Ptq (u). 

We will show later that in some cases the symmetrized operators will reduce to 
matricially free Gaussians. For instance, this is the case when the dimension matrix D 
is of dimension at least two and has one zero on the diagonal. When we evaluate the 
limit distribution of certain rectangular symmetric blocks under the associated partial 
trace, we obtain the case studied by Benaych- Georges [6]. 

In the case when U consists of one element, the corresponding array of symmetrized 
Gaussian operators was defined in [23] and thus we only have to extend that definition 
to the more general situation studied in this paper. 

Definition 3.1. If J is symmetric, then the symmetrized creation operators are opera- 
tors of the form 

$ {u) = { PM ^ P = Q 

for any (p, q) e J and u e U. Their adjoints p^ q (u) will be called symmetrized annihi- 
lation operators. The symmetrized Gaussian operators are sums of the form 

fipM = PpM + pp, q ( u ) 

for any (p,q) e 3 and u e U. Each matrix (u) Ptq (u)) will be called a symmetrized 
Gaussian matrix. 

We have shown in [23] that one array of symmetrized Gaussian operators gives the 
asymptotic joint distributions of symmetric random blocks of one Gaussian random 
matrix if the matrix D is non-singular. In other words, symmetric blocks of a square 
Hermitian Gaussian random matrix behave as symmetrized Gaussian operators when 
the size of the matrix goes to infinity. We will generalize this result to the case when 
we have a family of matrices indexed by a finite set U and we will allow the matrix D 
to be singular. The second situation leads to the following definition. 

Definition 3.2. If oo Ptq (u) ^ and oo qtP (u) ^ 0, then u) Ptq (u) will be called balanced. If 
UJ q,p( u ) = and 0J Pt q(u) ^ or u p<q (u) = and uu qiP (u) 0, then Q Ptq (u) will be called 
unbalanced. If uj pa {u) = u> qjP (u) = 0, then u) Ptq (u) = will be called trivial. 
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The main advantage of using symmetrized Gaussian operators is that in the case when 
the corresponding matrices B(u) are symmetric and D is non-singular, we can replace 
the free-boolean convolution of semicircle and Bernoulli laws by the free convolution of 
semicircle laws 

%(l)ffl%(2)ffl...ffl(T M (m) 
for p q (Lemma 6.2). The disadvantage is that this situation is less general and 
that is why it is the matricial freeness which remains the main notion of independence 
in the asymptotic regime. Although it is possible to define the notion of symmetric 
matricial freeness, it is derived from matricial freeness and is less universal. In fact, its 
formulation in [23] needs to be corrected in order to charaterize the algebra generated 
by symmetrized Gaussian operators and this modification also contributes to some 
loss of generality. The erratum to the definition of the symmetrically matricially free 
array of units is given in the Appendix. Nevertheless, the simplicity of Proposition 3.5 
and Corollary 3.1 and the importance of the underlying asymptotic model for random 
matrices seem to provide convincing arguments to treat this case. 

Example 3.1. In order to see the difference in computations of moments between the 
matricially free Gaussian operators and their symmetrizations, consider the following 
examples: 

*i« 2 ) = o 

= *l(pS,lP2,lP2,lP2,l) + *l(p|,lP*2Pl,2p2,l) = b ll + b h2h,l 
= *l(pl,2Pl,2Pl,2Pl,2) + *l(p*2p2,lP2,lPl, 2 ) = b l,2 + &1.2&2.1 

where, for simplicity, we took oo Ptq (u) = u PA for some fixed u in all formulas. If, for 
instance, a>i 2 = 0, then the third moment reduces to b\ x and the remaining three 
moments vanish. 

In order to state the relation between the symmetrized creation and annihilation 
operators, let us define the collective symmetrized units in terms of the array (l p , g ) of 
Definition 2.5 as 

respectively, for any (p, q) e J. The symmetrized units are described in more detail in 
the propositions given below. 

Proposition 3.1. Let J7 be symmetric and let Af PtQ be defined as in Definition 2.5 for 
any (p,q) e J. 

(1) If(q,q) e J, then l q , q is the orthogonal projection onto 

k 

(2) If (p, q) e 3 , where p q, then is the orthogonal projection onto 

cn p e cn g e 0(Ar p , fe e M q , k ). 
k 
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Proof. The above formulas follow directly from Definition 2.5. 



Proposition 3.2. The collective symmetrized units satisfy the following relations: 

( if {p,q}n{k,l} = 
l P ,?lfc,i = S if kandq = l 

{ %, q if (p, q) = (k, I) 
Proof. These relations are immediate consequences of Proposition 3.1. ■ 

Proposition 3.3. The following relations hold for any u,t e U: 

P P , q (u)p k ,i(t) = and p^ q (u)p* kl (t) = 
when {p, q] n {k, 1} = 0. Moreover, 

P* P , q (u)p k ,i(t) = 

when {p, q} ¥^ {k, 1} or u ^ t. Finally, if the matrix (b Piq (u)) is symmetric, then 

P*p, q (u)p P)q (u) = b Piq (u)l Piq 

for any (p, q) e J. 

Proof. These relations follow from Proposition 2.2. ■ 

In view of the last relation of Proposition 3.3, let us assume now that the matrices 
(bi,j(u)) are symmetric for all u e U, where 11 is finite. We would like to define a 
symmetrized array of algebras related to the array (Aij) of Proposition 2.3. For that 
purpose, we introduce collective symmetrized creation operators 

ueu 

and their adjoints p*^ q for any (p, q) e J . In turn, the sums 

Zp,q = Pp,, + P* g = ^ Up, q (u) 

ueu 

will be called the collective symmetrized Gaussian operators. 
Proposition 3.4. If (b Ptq (u)) is symmetric for any u e Zl, then 

Pl, q P P , q = K,<Ap,^ P P , q Pk,i = and pp tq p* ktl = 
whenever {p, q] n {k, 1} = 0, where b PA = T lu eu b p,q( u ) ■ 

Proof. The proof follows from Proposition 3.3. For instance, 

p*p,$pa = Ti pI,^ u )ppa^) = Ti b pM^i = b p,X, q 

u,teu ueu 

The remaining relations are proved in a similar way. ■ 

Let us discuss the notion of symmetric matricial freeness from an abstract point 
of view. We assume that the set J is symmetric, by which we mean that (j, i) e 
3 whenever e 3- Then we consider a symmetric array of subalgebras (Aij) 
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of a unital algebra A and we assume that (ljj) is the associated symmetric array of 
internal units. By I we denote the unital algebra generated by the internal units and 
we assume that it is commutative. Moreover, the array of states ((fiij) on A remains 
as in the definition of matricial freeness. Since the array (Aij) is symmetric, we thus 
associate two states, tpij and tpj ti , with each off-diagonal subalgebra Aij. 

Denote by J the set of non-ordered pairs {i, j}, even in the case when i = j. Instead 
of sets 3 m , we shall use their symmetric counterparts, namely subsets of the m-fold 
Cartesian product J x J7 x . . . x J7 of the form 

3m = {({il,3l}, {imJm}) ■ {h,Jl} # . . . # {im,jm}}, 

and their subsets 

%. m = {({iiJi}, • • • , {i m Jm}) ■ {k, jk} n {ifc+i, Jfc+i} ^ for 1 k < m - 1} 

where m e N. These sets comprise tuples of non-ordered pairs, in which neighboring 
pairs are different (as in the case of free products) and are related to each other as 
in non-trivial multiplication of symmetric blocks of matrices. It will be sometimes 
convenient to use the abbreviated notation zv k = {ik,jk}- 

In order to define the notion of symmetric matricial freeness, we shall first define the 
array of symmetrically matricially free units. This definition will differ from that given 
in [23] since we have discovered that the conditions on the units need to be strengthened 
in order that the symmetrized array of Gaussian algebras be symmetrically matricially 
free as claimed in [23, Proposition 4.2]. 

Definition 3.3. Let (Aij) be a symmetric array of subalgebras of A with a symmetric 
array of internal units (ljj), and let (<Pij) be an array of states on A defined by the 
family (<fj). If, for some a e A®j, where i j, it holds that 



or 



v(bl--a) = I V3iha) iU ^ 3 
WW) jo if i = j 

( PA° L *,* a ) jo if % # j 



for any b e A®j, then we will say that a is odd or even, respectively. The subspaces 



of A)j spanned by even and odd elements will be called even and odd, respectively. 



If each off-diagonal A® j is a vector space direct sum of an odd subspace and an even 
subspace, the array (Aij) will be called decomposable. 

Example 3.2. The idea of even and odd elements comes from the algebras generated 
by the symmetrized creation operators. If p ^ q and k is odd, it is easy to see that 

^ q(b^-p,pPp,q) = ^ q(^Pp,qPq,p ■ ■ ■ Pp,q) 
= * q (bPp,q) 

^Mq,qPpJ = 

for any b e A^ q and thus is odd. In turn, if k is even, 

%( b lq,q$,q) = ^ ' l( h Pq,pPp,q • • • Pp,q) 
= * q (bpp,q) 
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here is that only alternating products of these two operators do not vanish and thus a 
non-trivial product of odd or even order maps fl q into the linear span of simple tensors 
which begin with e PA (u), or e q ^ p (u), where u e 11, respectively. 

Definition 3.4. We say that the symmetric array (lij) is a symmetrically matricially 
free array of units associated with a symmetric decomposable array (Aij) and the array 
if for any state ipj it holds that 

(1) (pj(uiau 2 ) = (pj(ui)ipj(a)ipj(u2) for any as A and u±, u 2 e X, 

(2) if a r s A irjjr n Kery9 irjV , where 1 < r < m and (z^, . . . , %v m ) e % m , then 



for any as A, whenever one of the following cases holds: 

(a) w i n zif 2 iv 2 n z^ and a 2 is off-diagonal and odd, 

(b) ! n 2 = 2 n zf 3 and a 2 is diagonal or off-diagonal and even, 

where we set zv 3 = {j,j} for m = 2, and for any other (zfi, . . . , o/ m ) e J m the 
moment vanishes. 

Definition 3.5. We say that a symmetric decomposable array (Aij) is symmetrically 
matricially free with respect to (fij) if 

(1) for any a k e Ker(p iktjk n A fc ,j fc , where e [m] and (z^, . . . , z^ m ) e 3 m and for 
any state ipj it holds that 

ipj(aia 2 ...a m ) = 

(2) (ljj) is the associated symmetrically matricially free array of units. 

The array of variables (a it j) in a unital algebra A will be called symmetrically matri- 
cially free with respect to ((pij) if there exists a symmetrically matricially free array of 
units (ljj) in A such that the array of algebras (Aij), each generated by a it j + and 
ljj, respectively, is symmetrically matricially free with respect to (tfij)- The definition 
of *-symmetrically matricially free arrays of variables is similar to that of *-matricially 
free arrays. 

Proposition 3.5. Assume that 3 and (a Ptq (u)) are symmetric and a Ptq (u) > for any 
(p, q) e 3 and us U. Let 



and let (^ Ptq ) be the array of states on B(Ai) defined by the states . . . associated 
with the vacuum vectors Q 1 , . . . , Q r , respectively. Then 

(1) the array (A p , q ) is symmetrically matricially free with respect to (^! p ^ q ), where 
the associated array of units is 

(2) the \l> P)Q - distribution of u PA (u) is the semicircle law of radius 2a Ptq (u) for any 
(p, q) s J and us 11. 

Proof. The proof of (2) is the same as in the case when 11 consists of one element 
[23, Proposition 8.1], so we will only be concerned with proving (1). The proof of 
that one resembles that of Proposition 2.3. In that connection, let us observe that 
lp,<j ^ A Ptq for any (p, q) s J by Proposition 3.4, so that is why we did not include 



ip j (ali ltjl a 2 ■ ■ ■ a rn ) = tpj(aa 2 ■ ■ ■ a m ) 
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the symmetrized units in the set of generators. In order to prove symmetric matricial 
freeness of the array (A Pjq ), one needs to prove that the conditions of Definitions 3.4 
and 3.5 are satisfied. If p q, in view of commutation relations of Proposition 3.3, any 
element a e A® q := A PA n Ker^ p q is a linear combination of products of the form 

Pp,q(Ul) ■ ■ ■ Pp,q(Ui)vZ, q (tl) ■ ■ ■ Pp, q {tj), 

where ui, . . . , Ui, ti, . . . , tj e U and i + j > 0. Therefore, any element of A° pq of the 
above form, when acting onto a simple tensor w which does not begin with e Pjq (u) or 
e q,p( u ) f° r an y u i either gives zero or 

e Pim ( u i) ® e P2m (u 2 ) . . . ® e puqi (ui) ® w 

whenever w begins with e qu k{u) for some k and u, where all pairs (pj,qj) belong to 
{{Pi q)i (QiP)} an d alternate, where 1 < j < i. In turn, the result of the action of 
Ap q for p = q is identical to that treated in the proof of Proposition 2.3. Therefore, 
the action of the product a x . . . a m satisfying the assumptions of Definition 3.5 onto fl q 
either gives zero or a vector which is orthogonal to CQ q , which implies that 

#g(ai ...a m )=0, 

which proves the first condition of Definition 3.5. Let us prove now that (l p , q ) satisfies 
the conditions of Definition 3.4. The first condition holds since an analogous condition 
holds for the array (l p , g ). To prove the second condition, first observe that each off- 
diagonal Ap q is decomposable in the sense of Definition 3.3. Namely, the odd (even) 
subspace of A^ q is spanned by products of creation and annihilation operators of the 
form given above in which % is odd (even). Suppose first that 

{Pi,?i} n {P2,g2} = {P2,q2} n {p 3 ,q 3 }, 

with p 2 7^ q-ii for instance q\ = p 2 = ps, where indices Pk,qk are used instead of ik,jk- 
If a 2 e Ap 2 q2 is even, then a 2 . . . a m Q q is a linear combination of vectors of the form 
e P2tq2 (u) (x) w, where u e U and to is a simple tensor, which are left invariant by the 
unit l pi , qi by Proposition 3.1. If a 2 is odd, then a 2 . . .a m Q q is a linear combination of 
vectors of the form e q2<P2 (u) (B)w, where u e U and w is a simple tensor, which are killed 
by the unit l pim since q\ = p 2 ^ q 2 . Next, suppose that 

{Pi,9i} n {^2,^2} # {P2,Q2} n {P3,gs}, 

with p 2 q2, for instance q\ = p 2 and q 2 = p 3 . If a 2 is odd, then the vector a 2 . . . a m fl q 
is a linear combination of vectors of the form e P2m {u) (x)w, where u e 11 and to is a 
simple tensor, which are left invariant by the unit l pi , gi by Proposition 3.1. If a 2 is even, 
then a 2 . . . a m Q q is a linear combination of vectors of the form e q2tP2 (u)(S)w, where u s U 
and to is a simple tensor, which are killed by the unit l pim since qi = p 2 ¥= q 2 - Finally, 
if p 2 = q 2 and p 2 e {pi,qi}, then l Puqi leaves a 2 . . . a m fl q invariant. This completes the 
proof of symmetric matricial freeness. ■ 
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Corollary 3.1. Under the assumptions of Proposition 3.5, if J = [r] x [r], then the 
array of distributions of (u) M (u)) , u s 11, in the states (^ p , q ) is given by 

/ (u) cr li2 (u) . . . a lir (u) \ 

~, ^1,2 (u) a 2 , 2 (u) ... <T 2 , r (u) 

[a(u)\ = 

\ o- hr (u) a 2ir (u) . . . a r . r (u) J 
where o~ Pjq (u) is the semicircle law of radius 2a Pjq (u) for any (p, q) e 3 and us U. 

Each symmetric matrix of the above form plays the role of a symmetrized matricial 
semicircle law. It clearly differs from the matricial semicircle law of Corollary 2.1, but 
it contains all information needed to compute all mixed moments of the symmetrized 
Gaussian operators in the states ^i, . . . , \l/ r . For that reason, we use square matrices 
since we need two distributions of each off-diagonal operator u PA {u), in the states ^ p 
and ^ q . 

Remark 3.1. In the case when not all matrices A{u) = (a Ptq (u)) are symmetric, sym- 
metric matricial freeness is no longer sufficient due to more complicated relations be- 
tween the collective creation and annihilation operators, namely 

Pp,qPp,q = bp,q r p,q + "q,p r q,p 

where b p>q = Tj u eu b p,q( u ) and Kv = Tj u euKp( u ) for an y P ^ Q- Moreover, we have 
relations 

r p,qPp,q = Pq,pi Pp,q r p,q = Pp,qi Pp,qPp,q = ^p,q r p,qi lp,<? = r p,q r q,p 

and thus, assuming that b PA 7^ and b q ^ p ^ 0, we can see that the *-algebra A p>q gen- 
erated by p Ptq coincides with the *-algebra generated by p Ptq , p q , p ,r Ptq ,r qtP and thus it 
contains l PiQ . One can show that the mixed moments of interest, like those computed 
in Proposition 4.3, satisfy the same recurrences as those in the definition of symmetric 
matricial freeness. However, there are problems with mixed moments containing vari- 
ables like r Ptq and thus, on the level of algebras, we are not able to fit the array (A Ptq ) 
into the scheme of symmetric matricial freeness. This might indicate that from the 
point of view of independence the scheme of matricial freeness is the most natural one 
in this case. 

Let us only remark in this context that in this non-symmetric case the distributions of 
the operators u p<q (u) in the states ^ p<q can be easily determined. For that purpose, we 
shall need the probability measure i9 on K. corresponding to the two-periodic continued 
fraction (a, b, a, . . .). Its Cauchy transform is 

. z 2 + b - a - J{z 2 - b - a) 2 - Aab 

° {Z) = 2zb • 

It has the absolutely continuous part 

Znbx 

supported on \-\Ja — y/b\ ^ | x\ ^ yfa + Vb and an atom of mass 1/2 — a/2b at x = if 
a # b [21, Example 9.2]. In particular, if a = b, then is the semicircle law of radius 
2a. 
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Corollary 3.2. If J is symmetric and a p>q {u) > for any (p,q) e J7 = [r] x [r] and 

u e 11, then the array of distributions of (cu p<q (u)), u e 11, in the states (^ p>q ) is given 
by 

/ (Ti,! («) #1,2 («) ■ ■ ■ til,r( U ) \ 
~, v, #2,1 («) ^2,2(m) ••• #2,r( U ) 

[a(u)] = 

\ # r ,i(ti) # r , 2 («) • • • O r , r (w) / 

where and"& Ptq (u) is the distribution corresponding to the two-periodic continued fraction 
with the Jacobi sequence (b Ptq (u), b qtP (u), . . .) for any (p, q) e J and u e 11. 

Proof. The diagonal distributions follows from Corollary 2.1, whereas the off-diagonal 
ones will be computed in Example 5.2. ■ 

Each array [o(u)] is the symmetrized counterpart of [o(w)] and is called a sym- 
metrized matricial semicircle law. Both arrays uniquely determine the family of limit 
distributions of blocks of independent Hermitian Gaussian random matrices under the 
partial traces and thus under the trace. Moreover, there is a relation between their 
off-diagonal distributions of the form 

*W«) = «p,«(«) rn 

where □ is the s-free convolution defined by the subordination functions [21]. However, 
both arrays correspond to different matricial convolutions and require independent com- 
putations, although they are both based on matricial freeness in one way or another. 
In turn, the arrays (u Ptq (u)) and (Q Piq (u)) are different matricial generalizations of free 
semicircular families. Further, in the case when oo Ptq (u) = 0, the corresponding measures 
a Pt q(u), Kp,q{v) and $ P , q (u) should be replaced by Sq. 



4. Combinatorics of mixed moments 

It can be seen from the computations in Section 2 that the mixed moments of mate- 
rially free Gaussian operators are related to non-crossing pair partitions. The difference 
between them and free Gaussians is that one has to use colored non-crossing pair par- 
titions, where coloring is adapted to the pairs of matricial indices from J and to the 
additional indices from the set 11. It will be convenient to assume here that 11 = [i], 
where t is a natural number. 

For a given non-crossing pair partition n, we denote by B(n), C(n) and lZ(ir) the 
sets of its blocks, their left and right legs, respectively. If 7Tj = {l(i),r(i)} and iij = 
{Kj)i r U)} are blocks of 7r with left legs and and right legs r(i) and r(j), 
respectively, then 7Tj is inner with respect to tyj if I (J) < l{i) < r(i) < r(j). In that 
case iTj is outer with respect to 7Tj. It is the nearest outer block of 7Tj if there is no block 
n k = {l(k), r(k)} such that < l(k) < l(i) < r(i) < r(k) < r(j). Since the nearest 
outer block, if it exists, is unique, we can write in this case 

■Kj = ofc), l(j) = o(l(i)) and r(j) = o(r(i)). 

If 7Tj does not have an outer block, it is called a covering block. In that case we set 
o(iYi) = 7r , where we define 7r = {0,m + 1} and call the imaginary block. 
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Example 4.1. Let a e MC\ be as in Fig. 1. Its blocks are {1, 6}, {2, 3}, {4, 5} and 
the imaginary block is {0,7}. The left and right legs of tt are C{a) = {1,2,4} and 
7Z(cr) = {3,5,6}. The block {1,6} is the nearest outer block of both {2,3} and {4,5} 
and the imaginary block {0, 7} is the nearest outer block of {1,6}. 

Computations of mixed moments of (u Ptq (u)) in the states ty q are based on the classes 
of colored non-crossing pair partitions adapted to ordered tuples of indices. Since, 
in addition to the pair of matricial indices (p,q), we have an additional index u as 
compared with the case studied previously, we distinguish the matricial indices from 
the non- matricial one and we use an abbreviated notation for the pair v = (p, q). If the 
set 11 consists of one element, we recover the analogous definition of [23]. 

Definition 4.1. We will say that tt e MC 2 m is adaptedto the tuple ((z^, Ui), . . . , (v m , u rn )), 
where v k = (p k ,qk) e [ r ] x [ r ] an d u k g [t] for any k, if 

(a) (vi,Ui) = (vj,Uj) whenever {i, j} is a block of n, 

(b) qj = p (j) whenever {i, j} is a block of tt which has an outer block. 

The set of such partitions will be denoted by NC 2 m ((vi, u{), . . . , (v m , u m )). In turn, by 
NC 2 m q ((vi, ui), . . . , (v m ,u m )) we will denote its subset, for which q k = q whenever k 
belongs to a covering block. 

In order to find combinatorial formulas for the moments of our Gaussian operators, 
we need to use colored non-crossing pair partitons. It will suffice to color each tt g NC 2 m , 
where m is even, by numbers from the set [r]. We will denote by F r (n) the set of all 
mappings from the set of blocks of tt into [r] called colorings. By a colored non-crossing 
parition we will understand a pair (tt, f), where 7r e AfC m and / e F r (n). The set 

B(*,f) = {(*!, /),..., (7T fc ,/)} 

will denote the set of its blocks. We will always assume that also the imaginary block 
is also colored by a number from the set [r] and thus we can speak of a coloring of n. 

Example 4.2. If it g J\fC 2 mq ({v u m), . . . , (v m , u m )), then ((^, m), . . . , (v m , u m )) defines 
a unique coloring of n and tt in which the block containing k is colored by p k for any 
k and the imaginary block is colored by q. Consider the first colored partition given in 
Fig. 1, where, for simplicity, we assume that t = 1 and thus we can skip it's. If we are 
given the tuple of pairs 

(^2, z^ 4 ) = ((2, 1), (2, 2), (2, 2), (2, 1)) 

then tt is adapted to this tuple. If q = 1, then the unique coloring of tt defined by the 
given tuple and the number q = 1 is given by }\ since the imaginary block must be 
colored by 1 and the colors of blocks {1,4} and {2,3} are obtained from the numbers 
assigned to their left legs, i.e. /i(l) = 2 and /i(2) = 2. 

Definition 4.2. Let a real-valued matrix B{u) = (b it j(u)) e M r (R) be given for any 
us [t]. We define a family of real- valued functions b q , where q e [r], on the set of 
colored non-crossing pair-partitions by 

b q (TT, f) = &,(7Ti, /) . . . b q (n k , f) 

where tt g NC 2 m and / e F r {rr), and where b q is defined on the set of blocks B(tt, /) as 

b q (TT k ,f) = bij(u), 
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FIGURE 1. Colored non-crossing partitions 

whenever block ix k = {r, s} is colored by i, its nearest outer block is colored by j and 
u r = u s = u, where we assume that the imaginary block is colored by q. 

It should be remarked that in this paper b q (7r, f) may be equal to zero even if tt g 
A/"C^((z>i, Ui), . . . , (v m , u m )) since we assume that the matrices B[u) may contain zeros. 
Let us also recall our convention saying that if MC 2 m or its subset is empty we shall 
understand that the summation over tt g A/"C^ or over this subset gives zero. In 
particular, this will always be the case when m is odd. 

Proposition 4.1. For any tuple ((vi, ui), . . . , (v m , u m )) and q g [r], m s N, where 
Vk = (Pk, Qk) e [r] x [r] and g [t] for each k, it holds that 

* ? K li91 (« 1 )...w Mm (u m ))= J] b q (n,f) 

TTeMCl 2tq ((vL,U 1 ),...,(vm,Um)) 

where f is the coloring of 7r defined by ((tfi, ui), . . . , (v m , u m )) and the number q. 

Proof. The proof is similar to that of [23, Lemma 5.2] and it reduces to showing that 
if ei, . . . , e m g {1, *} and tt g HC 2 m q {{vi, U\), . . . , (v m , u m )), where m is even, then 

VM\,J U J ■ ■ ■ P e p2, q J u ™)) = b q (n, f) 

where / is the coloring of tt defined by the collection of indices {pk, k g C(ty)} associated 
with the left legs of the blocks of tt and the index q coloring the imaginary block. The 
only difference between the former proof and this one is that to each block of (tt, f) we 
assign a matrix element of B(u) for suitable u (u is the same for both legs of each block 
since the partition satisfies condition (a) of Definition 4.1). ■ 

Using Proposition 4.1, we can derive nice combinatorial formulas for the moments of 
sums of collective Gaussian operators 

where 

ueu 

and for that purpose, denote by NC^ q [r] the set of all non-crossing pair partitions of 
[m] colored by the set [r], with the imaginary block colored by q. 
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Lemma 4.1. The moments of u in the state ty q , where q e [r], are given by 
where b q (ir, f) = b q fa, f) . . . b q (n s , f) for n = fa, . . . , iv s } and 

bqfaj) =YjKj( U ) 
u 

whenever block n k is colored by i and its nearest outer block is colored by j. 
Proof. Using Proposition 4.1, we can express each summand in the formula 
*«(w m ) = 2 ^K ll5l (Mi)--.w Mm W) 

Pi,qi,ui,...,p m ,q m ,u m 

in terms of b q (n,f), where 7r e MC^ q {{vi, ui), . . . , (z/ m , u m )), with = (p fc ,g fe ) and 
/ is the coloring defined by the tuple ((vi, ui), . . . , {v m , u m )). A different choice of 
((z>i, Mi), . . . , (v m , u m )) must lead to a different collection of colored partitions (ir,f) 
since even if the same partition appears on the RHS of the formula of Proposition 
4.1, the coloring must be different. In fact, 7r is uniquely determined by the sequence 
e = (ei, . . . , e rn ) which appears in the nonvanishing moment of type 

^q(p £ p i qi fa)... P ;Z,qJ U m)), 

where e k e {1, *}. Moments of this type are basic constituents of ^ q {u m ). For instance, 
ife=(*,*,l,*,l,l), the associated unique non-crossing pair partition is a of Fig. 1. If we 
keep 7r and change at least one index in the given tuple to which n is adapted, we either 
get a tuple to which n is not adapted (if we change q k for some k e since this index 
is determined by the color of its nearest outer block) or we obtain a different coloring / 
(if we change p k for some k e Cfa). Therefore, when we add all contributions b q (7r, f) 
of the form given by Proposition 4.1 corresponding to different tuples of indices, we 
will obtain contributions corresponding to different elements of MC 2 m [r]. Therefore, 
in order to prove the formula for ^ q (u m ), we only need to justify that all colored 
partitions from MC 2 m q [r\ really do contribute to ^ q (u m ). Thus, let (n,f) e MC 2 m q \r\ 
be given. There exists a unique e = (ei, . . . , e m ) associated with ir. In turn, / determines 
((pi, qi), . . . , (p m , q m )) by an inductive procedure with respect to blocks' depth. Thus, 
if {k,l} is a covering block, we choose q k = qi = q and p k = pi = f(k) = /(/). Next, 
if {k,l} = o({i,j}), then we choose qi = qj = p k and pi = pj = f{i) = f(j), etc. We 
proceed in this fashion until we choose all p k , q k , k e [m\. This completes the proof. ■ 

Example 4.3. For square arrays (ujij(u)) of dimension r, where u e U, and for any 
q s [r], we have 

= J] ^MA^^Plq^PUt)) + * 9 (^ 9 (s)P^-(*)Pij(*)Pi,«(s))) 
i,j,s,t 

i,j,s,t 

^ b q faf)+ J] b q (a,f) 

feF r (n) geF r (a) 
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where tt = {{1, 2}, {3, 4}} and a = {{1, 4}, {2, 3}} and where b q (ir,f) and b q (a, /) are 
defined as in the proof of Lemma 4.1. 

Note that u coincides with the sum of all symmetrized Gaussian operators 

and thus Lemma 4.1 gives a formula for the moments of u as well. However, this is 
not the case for the mixed moments of symmetrized Gaussian operators (u Ptq (u)) in the 
state ty q . These are based on the class of colored non-crossing pair partitions adapted 
to ordered tuples of indices of type (tv , u), where w is an abbreviated notation for the 
set {p, q}. 

The definition of this class is given below. Note that it is slightly stronger than 
that in [23], which is related to the stronger definition of symmetric matricial freeness 
discussed in Section 3. 

Definition 4.3. We say that tt e MC 2 m is adapted to the tuple ((u>i, Ui), . . . , (rv m , u, m )), 
where rvk = {pk,Qk} and (pk,qk,Uk) e [r] x [r] x [t] for any k, if there exists a tuple 
((v 1 , iti), . . . , (v m , u m )), where v k e {(p k , q k ), (qk,Pk)} for any k, to which tt is adapted. 
The set of such partitions will be denoted by MC m ((wi, Ui), . . . , (zv m , u m )). Its subset 
consisting of those partitions for which n s J\fC^ q ((vi, Ui), . . . , (v m , u m )) will be denoted 

Example 4.4. For simplicity, consider the case when t = 1 and we can omit it's. If we 
are given the tuple of sets 

«*,«/ 4 ) = ({1,2},{1,2},{1,2},{1,2}), 

then the partition 7r of Fig.l is adapted to it since there are two tuples, 

((2,1), (1,2), (1,2), (2,1)) and ((1, 2), (2, 1), (2, 1), (1, 2)), 

to which 7T is adapted. If q = 1 or q = 2, then the associated coloring is given by / 2 or 
/ 3 , respectively. Thus, there are two colored partitions, (tt, f 2 ) and (tt, / 3 ), associated 
with the given partition 7r and the tuple (zfi, a^, a^, w±). Nevertheless, once we choose 
the coloring of the imaginary block, the coloring of tt is uniquely determined. In turn, 
if we are given the tuple of sets 

^4, = ({1, 2}, {1, 2}, {1, 2}, {2, 2}, {2, 2}, {1, 2}) 

then a of Fig.l is adapted to it since it is adapted to the tuple 

((2,1), (1,2), (1,2), (2,2), (2,2), (2,1)). 

Here, there is no other tuple of this type to which o would be adapted and therefore the 
only coloring associated with tt is given by g and the only color which can be assigned 
to the imaginary block is 1. 

Example 4.5. Let us give an example of a partition which satisfies the conditions of 
Definition 2.3 of [23], but is not adapted in the sense of Definition 4.3 of this paper. 
Namely, the partition 



tt = {{1, 8}, {2, 7}, {3, 6}, {4, 5}} with diagram 











n 
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is not adapted to ({1, 1}, {1, 2}, {1, 2}, {2, 2}, {2, 2}, {1, 2}, {1, 2}, {1, 1}) since there is no 
tuple of ordered pairs (v±, . . . , v%) to which 7r would be adapted according to Definition 
4.1 ({1, 8} must be colored by 1, this forces {2, 7} to be colored by 2 and thus {3, 6} to 
be colored by 1 by condition (b) of Definition 4.1, but 1 ^ {2, 2}). At the same time, n 
satisfies the conditions of Definition 2.3 of [23]. On the other hand, the partition 



a = {{1,6}, {2, 5}, {3, 4}} with diagram 

is adapted to ({1, 1}, {1, 2}, {2, 2}, {2, 2}, {1, 2}, {1, 1}). The actual reason why a is 
adapted to the associated tuple (and 7r is not) is that the pair {1,2} appears in this 
tuple an odd (even) number of times in a row and the associated legs belong to different 
blocks. This feature is related to Definition 3.4, in which even and odd elements ap- 
pear, and is a consequence of the difference between the action of odd and even positive 
powers of symmetric off-diagonal blocks T p q of a matrix onto basis vectors. Namely 

W^p and T p k q (V q )^V q 

for j odd and k even, respectively, where V r = span{ej : i e N r } and the dependence of 
all objects involved on n and u is supressed in the notation. 

It can be seen from Example 4.4 that if it e AfC^^wi, Ui), . . . , (w m , u m )), then there 
may be more than one colorings of n defined by the colorings of the associated tuples 
{{vy, ui), . . . , {v m , u m )). Therefore, these tuples may produce more than one coloring 
of 7r defined by the sets of matricial indices. However, when we fix q and require that 
7r e J\fCm >g ((vi, Ui) ■ ■ ■ , (vm, u m)), we obtain a unique coloring of rr since we have a 
unique associated tuple ((vi,Ui),. . ., (v m ,u m )). 

Proposition 4.2. If ir <= MC^^zVi, Mi), . . . , (rv m , u m )), there is only one associated 
tuple ((vi^Ui), (v m ,u m )) for which n e MC 2 ni:q ((v u u x ), . . . , {v m ,u m )). 

Proof. If {m — 1, m} is a block, then the second index of v m -\ = v m must be q and 
the imaginary block is also the nearest outer block of the block containing m — 2. This 
allows us to treat the partition n' obtained from 7r by removing the block {m — 1, to} in 
the same way, which gives the inductive step for this (same depth) case. If {m — 1, m) 
is not a block, then the block containing m is the nearest outer block of that containing 
m — 1 and thus the second index of v m -i must be equal to the first index of v m and 
the second index of v m is q. This determines v m and v m -i and gives the first inductive 
step for this (growing depth) case. Proceeding in this way, we determine vi, . . . , v m in 
a unique way. ■ 

Proposition 4.3. For any tuple ((zifi, Ui), . . . , {w m , u m )) and q e [r], where u>k = 
{Pk, Qk] o,nd pk,qk e [r], Uk e [t] for any k and m e N, it holds that 

^(^ ll(!1 («i)---^ m , gm W)= ^ b g (ir,f) 

T6A/"C^ ig ((a^,ui),...,(7f m ,w m )) 

where f is the coloring of n defined by ((u>i, Ui), . . . , (zv m ,u m )) and by the number q. 

Proof The LHS is a sum of mixed moments of the type computed in Proposi- 
tion 4.1 with perhaps some p^s interchanged with the corresponding q^s. These mo- 
ments are associated with tuples of the form ((i> 1; ui), . . . , {v m , u m )). It follows from the 



n 
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proof of Proposition 4.1 that with each moment of that type we can associate the set 
NC 2 rn q ((vi,ui), . . . , (v m , u m )) and each such moment contributes 

2 b q (*,f), 

1tE A f C 1 ,((n,''l),-,(%.»m)) 

where J\fC m>q ({vi,Ui), (vm,u m )) c A/"C^ g ((^i, m), (u/ m ,u m )). It is clear from 
Definition 4.3 that 

U NC m , q ((vi, Ml), • • • , («*n, Tim)) = AfC^, ((«* , Mi), ... , («/ m , Um)), 
(vi,...,v m ) 

where the union is taken over the tuples . . . , v m ) which are related to {rvi, . . . , zv m ), 
i.e. if u>k = {Pk,Qk}, then Vk e {(pfc,?fc), (?fc,Pfc)}- Since the sets on the LHS are disjoint 
by Proposition 4.2, the proof is completed. ■ 

Finally, we would like to introduce a subset of NC m q ((rOi, ui), . . . , (w m , u m )) consist- 
ing of those Ti for which the coloring / satisfies certain additional conditions. Namely, 
we would like to distinguish only those colored partitions (tt, f) whose blocks are not 
colored by q. Thus, only the imaginary block is colored by q, which is not a contra- 
diction since the imaginary block is not a block of n. This subset will be denoted by 
NCT mq ((w u Mi), ... , (w m , Um)). 

Example 4.6. In Fig.l, (tt, /i) e A r CX 4 2 1 ({l, 2}, {2, 2}, {2, 2}, {1, 2}) since q = 1 and no 
blocks of 7r are colored by 1. In turn, the remaining colored partitions are not of this 
type since the colors of their imaginary blocks ars assigned to other blocks, too. 

Proposition 4.4. For given q e [r], if the only unbalanced operators in the given 
moment are of the form Piq (u) = ui q . p {u) = oo Ptq (u), where q ^ p e [r] and u e [t], then 

M^WiM ■■■u Pm , qm {u m )) = Yi h q{^J) 

■KeMCX^ lq {(m,u-i),...,{ifm,u m )) 

where f is the coloring of n defined by ((z^l, ui), . . . , (zv mi u m )) and the number q. 

Proof. The difference between the considered mixed moment and that of Proposition 
4.3 is that all operators are balanced except those involving the index q. Namely, if 
q e {pk,qk} for some k, then we have u Pktq (uk) or ui qkA {uk) instead of the symmetrized 
Gaussian operator u Pkjqk (uk) for such k and Uk e [t\. Since the first index is these opera- 
tors is different from q, we have to eliminate from the set J\fC m ((zvi, u±) . . . , (rv m , u m )) 
those partitions in which q colors any blocks of tx since it is always the first index 
which colors the block. This means that we are left only with the contributions from 
7T g MCX mq {{zVi, Mi) ... , (w rn , iij), which completes the proof. ■ 

The formula of Lemma 4.1 holds irrespective of the number of trivial operators in 
the arrays (ou Piq (u)) except that some of the contributions vanish. We can derive a 
similar formula in the case when the only unbalanced operators are of the type studied 
in Proposition 4.4. By J\fCX m q \r\ we shall denote the subset of MC m q \r\ consisting of 
those partitions in which no blocks (other than the imaginary block, which formally 
does not belong to the partition) are colored by q. 
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Lemma 4.2. If the only unbalanced operators in the arrays (u)ij(u)) are of the form 
&p,q( u ) = £q,p( u ) = U} p,q( u )> where p 7^ q e [r] and ue [£], i/iera 

(^,/)eArcX2 i9 [r] 

where b q (ir, f) = b q (n u f) . . . b q (ir s , f) for n = {ty 1 , ...,tt s } and q e [r], meN, and 

b q (lTk,f) =Yi b i,j( U ) 
u 

whenever block ir k is colored by i and its nearest outer block is colored by j . 

Proof. The proof is similar to that of Lemma 4.1 (Proposition 4.4 is used). ■ 

Example 4.7. Consider square arrays (uij(u)) of dimension r, where u e U, which 
have unbalanced operators of the form 0Ji, q {u) = O0i, q (u) for any i ¥^ q, where q e [r] is 
fixed and us [t]. We have 

Q = 2 2 ^'( u ) = 2 2 

i¥=q 

and thus 

s,t,i,j 



2 + mow*)) 



s,t,i,j 



= ^ 6 9 (tt,/) + 2 

feF r , q {n) feF r , q (a) 

where F r ,s(C) denotes the set of colorings of £ by the set [r]\{g} and the remaining 
notations are as in Example 4.3. Note that some of the contributions of Example 4.3 
disappear. 

5. Hermitian Gaussian random matrices 

The context for the study of random matrices originated by Voiculescu [32] is the 
following. Let \x be a probability measure on some measurable space without atoms and 
let L = ni^ P <oo ^(aO be endowed with the state expectation E given by integration 
with respect to \x. The *-algebra ofnxn random matrices is M n (L) = L® M n (C) with 
the state 

r(n) = E (x) tr(n) 

where tr(n) is the normalized trace over the set {ej : j e [n]} of basis vectors of C n . 

In order to study the asymptotics of symmetric blocks of random matrices, we par- 
tition each set [n] into disjoint non-empty intervals 

I AM 

[n] = Ni u . . . u N r , where lim — > dj ^ 

n^cc n 

for any j, and we set again D = diag(di, c^, • • • > d r ) to be the associated diagonal 
dimension matrix. The numbers d±, . . . ,d r will be called asymptotic dimensions. 
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Definition 5.1. By a Hermitian Gaussian random matrix we shall understand a Her- 
mitian n-dimensional matrix Y{n) = {Y id (n)) of complex- valued random variables such 
that 

(1) the real- valued variables 

{ReY id (n), lmY^(n) : 1 ^ i < j < n) 

are independent Gaussian random variables, 

(2) E(Yij(n)) = for any i,j,n, 

(3) E(|y^(n)| 2 ) = v p Jn for any (i,j) e N p x N q and n, where U : = (v p , q ) e M r (R) . 

The asymptotic joint distributions of symmetric blocks of one Hermitian Gaussian 
random matrix in the case when all asymptotic dimensions are positive have been 
studied in [23], where we have shown that they have the same limit joint distributions 
as the symmetrized Gaussian operators. In this paper, we shall generalize this result 
to the case when we have a whole family of Hermitian Gaussian random matrices 

{Y(u, n),u e 11} 

for any neN. Moreover, we will allow some asymptotic dimensions to vanish. 

Definition 5.2. We will say that the family of Hermitian Gaussian random matrices 
of Definition 5.1 is independent if the variables 

{ReY it j(u, n), lmYij(u, n) : 1 < % < j < n and u e 11} 

are independent and 

E(|y M ( M ,n)| 2 ) = ^M 
n 

for any e N p x N q and n e N,u e W, i.e. the variances of the complex- valued 
variables Y it j(u, n) are identical within each block defined by N p x N q for any ue 11. 

As in the case of one matrix, we will study the asymptotic joint distributions of the 
submatrices of the form 

T Ptq (u,n)= Y itj (u,n)®e itj (n), 

where {eij(n)\l ^ i,j ^ n} is a system of matrix units for each n and 

N p , q = (N p x iV ? ) u (N q x JV p ) 

for p, q e [r] and u e 11, called independent symmetric random blocks. Non-Gaussian 
random matrices [12] and the corresponding symmetric blocks can be treated in a similar 
way, but our proofs rely on the Gaussian case treated in [23,31]. 

We would like to find a Hilbert space realization of the limit joint distributions of 
independent symmetric random blocks under r(n) and under 

r q (n) = E (x) tr q (n), 

where tr q (n) is the normalized partial trace over the set {ej : j e N q }, where by normal- 
ization we understand division by |iV g |. 

Definition 5.3. The sequence of symmetric random blocks (T Pjq (u, n)) nG N, where (p, q) e 
J, u e 11 are fixed, will be called balanced if d p > and d q > 0, unbalanced if 
(d p = and d q > 0) or (d p > and d q = 0) and evanescent if d p = and <i g = 
(cf. Definition 3.2). 
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Example 5.1. Consider the sequence of Hermitian Gaussian random matrices of the 
block form 

Y(„\ - ( Sl M \ 

Y {U) ~ { St, 2 (n) S 2 , 2 (n) )> 

where both diagonal blocks are Hermitian for any n e N (u is omitted for simplic- 
ity). Thus, we have three symmetric blocks, two diagonal ones, Ti ;1 (n) = £1,1(71) and 
^2,2(71) = S 2j2 (n), and the off-diagonal one that can be identified with the Hermitian 
matrix 

S 1 , 2 (n) \ 
S* h2 (n) J ■ 

If we suppose that 

liVil \N 2 \ 
lim = d\ = and lim = d 2 > 0, 

then (T ltl (n)) ne ^ is evanescent, (T 2t2 (n)) ne ^ is balanced and (T lj2 (n)) ne ^ is unbalanced. 
It will follow from Theorem 5.1 that the moments of such symmetric blocks under 
T q (n) tend to the moments of trivial, balanced and unbalanced symmetrized Gaussian 
operators, cDi,! = 0, ^2,2 = ^2,2 and Qi t2 = ou 2t i, respectively, under ty q , where q e {1, 2}. 

In fact, it is easy to predict that if D is singular, then all mixed moments involving 
evanescent blocks will tend to zero. However, if we compute mixed moments involving 
unbalanced blocks consisting of pairs of mutually adjoint rectangular blocks whose one 
dimension grows too slowly, it is less obvious that one can still obtain nonvanishing 
limit moments under a suitable partial trace. 

In the Fock space realization of limit joint distributions of blocks under r(n) we will 
always use the convex linear combination of states . . . , \P r of the form 

r 

* = 2 d q %, 

q=l 

which is a state on B(A4), but in the study of the limit joint distributions of blocks 
under r q (n) we will use the single state ty q . In this case, it will not be convenient 
to assume that is a special case of ^ with d q = 1 and the remaining asymptotic 
dimensions vanishing. 

As we have proved in [23], the asymptotic joint distributions of blocks of one random 
matrix are expressed in terms of the matrix B = DV being the product of the dimension 
matrix and the symmetric matrix of numbers V = (v Ptq ) related to the variances, namely 

b p , q = d p v p , q , where — = E(\Yij(n)\ 2 ) 

for each complex variable Yij(n) in the block T Ptq (n) of the matrix Y{n). In the present 
generalization, when we deal with a family of independent random matrices, the situa- 
tion is similar and we introduce the family of matrices 

B{u) = DV(u), 

where V(u) = {v ps {u)) for any u e 11 and thus b Ptq (u) = d p v p>q {u) for any p,q e [r] 
and u e U. Let us observe that even in the case when the random matrices Y(u, n) are 
Hermitian, which implies that the matrices V(u) are symmetric, the associated matrices 
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B{u) are symmetric only if all dimensions are equal. In turn, if all variances are equal 
to one, then 





( d x 


d\ . 


■ d 1 \ 


B(u) = 




d 2 . 


. d 2 










V d r 


d r . 


. d r ) 



for any u e 11. Even this simple situation is interesting when we investigate joint 
distributions of blocks. Of course, when computing asymptotic distributions of matrices 
Y(n,u), we do not need to split them into blocks and free probability suffices. 

Theorem 5.1. Let [Y(u,n) : u e U) be a family of independent Hermitian Gaussian 
random matrices for any n s N and let (T p>q (u, n)), where u e U and n e N, be arrays 
of their symmetric blocks. Then 

lim T q (n)(T Puqi (u u n) . . . T Pmtqm (u m , n)) = * g (a5 Pl , gi (ui) . . . uJ Pm , gm (u m )) 

n—>cc 

for any (f>i, ft), . . . , (p m , q m ) e J, q e [r], and u ± , . . .u m e U, where u p , q {u) 's are sym- 
metrized Gaussian operators associated with matrices B{u). 

Proof. We would like to generalize the proof of [23, Theorem 9.1]. First, consider the 
partial trace T q (n)(T puqi (ui,n) . . .T Pmtqm (u m ,n)) for any q e [r] and even m = 2s. We 
need to enumerate contributions from moments of the form 

which are associated with 7r(7) e J\fC^ q ((zv i, Ui), . . . , (zv m , u rn )), where Wi = {pi,qi} 
for i e [m], defined by permutations 7 : [m] — > [m] such that 

7 2 = id, 7(fc) ^ k and L l{k) = j k , j l{k) = i k , u l{k) = u k 

Each such moment is 0{n~ s ~ r ) as n — > 00. For convenience, we restrict our attention 
to the case when 

mi jrn) ^ ^p m ^ ^9m' 

since the remaining cases can be treated in a similar way, with some p^s interchanged 
with the corresponding q^s. Moreover, we must have q = P\ since the action of 
Y imtil (u m ,n) onto vectors ej, where i\ $ N q and j e N q , will give zero. Assuming 
therefore that q = pi, we can see that the number of non-zero moments of the consid- 
ered type is 2 7 ©71(7), where 7 runs over the set of permutations of [to] such that 

7 2 = id, j(k) ^ k and q lik) = p k , p l{k) = q k , u^ k) = u k 

for any k e [to]. At this stage of the proof, let us observe that the main difference 
between this proof and that in [23] is that the set of permutations of [to] just described 
may be smaller since these must be adapted to u±, . . . ,u m , which, in general, leads to 
a smaller class of partitions, namely 

TV 1 , . . . , H) m ) 

for any U\, . . . ,u m . Here, one has to take into account the erratum given in the Ap- 
pendix which results from the stronger definition of partitions adapted to the tuples 
y Ifj, . . . , TV m ) (Definition 4.3). Thus, only if 71(7) e MC 2 ((wx, u-y), . . . , (zv m , u m )), 
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the number of independent indices which enter in the computation of 6^(7) is s + 1. 
However, the enumeration of the cardinalities of O n (7) for such 7 remains the same: 

n (7) = card{(ii, i rW , v (s) ) e N P1 x N Pr(1) x . . . x N Pr(s) }, 

where 

ft(7r( 7 )) = {r(l),r(2),...,r(s)} 

is the set of right legs of 77(7). Thus the dimension formula for O n (7) is again of the 
form 



©nw) n Pi n p r ii) ■ ■ ■ 



where n p = \N P \ for any p. For those independent indices which are pairwise different 
and correspond to partitions 77(7) e J\fC mq ((wi, ui), . . . , (w m , u m )), the corresponding 
expectation of complex Gaussian random variables K(Yi 1 j l (ui,n) . . ■ Yi m ,j m (u m ,n)) is 
the product of variances indexed by matricial indices and labelled by elements of 1Z, 
namely 

E(y JfcJfc K,n)y l7(fc)j7(fc) ( w ,n 7(fc) )) = v j^A 

where j(k) e TZ(tt('j)), with z 7 (&) = j^, j^^) = ik and u^ha = u^, and where the 
assumption that variances are equal within blocks is used. Let us conclude that in the 
computations of mixed moments of symmetric blocks taken from a family of independent 
matrices rather than from one matrix, we just have to replace the class of parititions 
MC 2 nhq (iOi, . . . , tu m ) by AfC^ ui), ... , («/ m , Um)) and the variances v PA by v p , q (u). 

We now collect expectations corresponding to all n e J\fC mq ((zt/i,ui), . . . , (zi> m , u m )) 
and take into account that 

hm = d p . ..dp 

n^cc n q n s w w 

where we use the fact that p\ = q and thus n Pl /n q = 1. Suppose first that all symmetric 
blocks in the considered moment are balanced. Then d k > for all relevant fc's, which 
implies that the contribution from each expectation corresponding to the given it is 
b q (ir, f), where / is the coloring of it defined by q and the tuple (zui, . . . , iv m ). Of course, 
the values of b q {n, f) depend on the matrices B(ui), . . . , B{u m ). Thus, we obtain 

2 b g (n,f) 

■KeJ\TC'i n>q ({vi,u{),...,(Vm,Um)) 

as n — > 00. Since this expression is equal to 

^(^WlK) ■ ■■^Pm,qm( U m)) 

by Proposition 4.3, our assertion is proved. In turn, if there is an unbalanced block, 
then the product of asymptotic dimensions given above is non-zero if the associated 
vanishing dimension does not label any inner blocks of it. This is possible if and only 
if 7T 6 A/'CI^ 9 ((n'i, Ui), . . . , («/ m , u m )). In this case, we obtain 

2 b q (*,f), 

),-.;{w m ,Um)), 

which, by Proposition 4.4, is equal to ^ q (u} puqi (ui) . . ■ w Pmiqjn (u m )) . Finally, if the in- 
vestigated moment contains one evanescent block, then in each product of asymptotic 
dimensions cL eL , , at least one dimension vanishes and thus the limit moment is 
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zero, which agrees with the corresponding moment of symmetrized Gaussian operators 
since one of them is trivial. This completes the proof. ■ 

Corollary 5.1. Under the assumptions of Theorem 5.1, 

lim T(n)(T puqi ( Ul ,n) . . . T Pmit]m (u m , n)) = *(tDp 1)gi (ui) . . . u Pm , qm {u m )) 

n—>oo 

for any (p 1 ,q 1 ),. . ., (p m , q m ) e J7 and U!,...u m e U. 
Proof. We have 

where \N q \/n — » d q as n — » oo for any q s [r], and thus the assertion follows from 
Theorem 5.1. ■ 

Example 5.2. Let (T p ^ q (n)) ne ^ be a sequence of off-diagonal balanced symmetric blocks 
for fixed (p, q) e J (u e 1/ is omitted for simplicity). By Theorem 5.1, 

lim r 9 (n)(T-(n)) = ), 

where we have 

ei,...,e m e{l,*} 

and only the summands which correspond to tuples (e±, . . . ,e m ) which define non- 
crossing pair partitions give non-zero contributions (cf. [23, Lemma 4.1]), where we 
skip u in all notations. Moreover, if such a tuple defines a non-crossing pair partition, 
then all blocks of depth one contribute b Ptq since the associated creation operators act 
on fl q , the blocks of depth two contribute b qjP since the associated creation operators 
act on e p>q , blocks of depth three contribute again b PA since the associated creation 
operators act on e qtP (x) e Piq , etc. By Proposition 4.3, for even m, we have 

q, (Qm) = y b \B {n)\ b \B e {n)\ 
*q\ M p,q) /j P.9 9.P 

since the set of non-crossing pair partitions of [m] which are adapted to the tuple 
{{P) ?})•••) {Pi q}) °f lenght m and to the number q can be identified with the set of 
all non-crossing pair partitions of [m], where B (n) and B e (ir) denote the blocks of ir 
of odd and even depths, respectively. Therefore, the Cauchy transform of the limit 
distribution can be represented as the two-periodic continued fraction with alternating 
Jacobi coefficients (b p>q , b q>p , b p>q , . . .). 

Example 5.3. Let (T p ^(n)) ne ^ be a sequence of off-diagonal unbalanced symmetric 
blocks, where p q are fixed and d q = and d p ^ (u e 11 is omitted for simplicity). 
If m is even, then 

Jim T q (n)(T™(n)) = %(^ q ) = %((p* p>q p P , q r /2 ) = 6^ 

as n — > oo since u qtP = 0. Of course, if m is odd, we get zero. The remaining partial 
traces vanish, in particular 

kT,( B )(^( B ))=*,(^)=0 

as n — > oo, since u Pjq Q p = for p ^ q. 
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It should be clear that if all variances are equal to one, all mixed moments of Theorem 
5.1 become polynomials in asymptotic dimensions. This case is interesting especially 
if not all asymptotic dimensions are equal. Although finding explicit forms of such 
polynomials may be highly non-trivial, a combinatorial formula can be given. 

Definition 5.4. For any tuple ((wi, iti), . . . , {u> m , u m )), where u>\ u . . . u w m c [r] and 
Ui, . . . , u m e U, with m even, and for any q s [r], let 

Q s (d u d 2 ,...,d r )= ^ n dj^w 

7T£A/'CS l ((jW,Ul),...,(o* n ,U m )) jewiu...uw m 

be a polynomial in asymptotic dimensions d 1 , . . . , d r , where TZjiir) is the set of right 
legs of tv which are colored by j. 

Remark 5.1. It is not hard to see that all these polynomials are homogenous and can 
be written in the form 

Q s (d 1 ,d 2 ,...,d r ) = ^ M ku ..., kr d^...d r k ' 

0^k 1 ,...,k r t : s 
ki + ... + k r — s 

where Mk lt ... t k r is the number of those it e MC 2 mq {fyW i, u\), . . . , (zv m ,u m )) which have 
kj right legs colored by j, where j e [r]. Clearly, the coefficients Mfc li ... i A; r depend on 
the tuple ((u>i,ui), . . . , {xv r , u r )) and on q and may vanish. 

Corollary 5.2. Under the assumptions of Theorem 5.1, if v PA {u) = 1 for any p,q e 
[r],u£ U, then for any m = 2s 

lim rg(n)(T pi> g 1 (ui,n)...Tp mj , m (ti m ,n)) = Q s {d l ,d 2 , . . . ,d r ) 

n— >co 

where Q s (d\, d 2 , ■ ■ ■ , d r ) is the polynomial of Definition 5.4- 

Proof. Non- vanishing contributions correspond to n 6 MC? m g ((^i, Ui), ■ ■ ■ ,(w m , u m )) 
and are equal to products of asymptotic dimensions, namely 

where TZ(tt) = (r(l), . . . ,r(s)}, which proves our assertion. ■ 

In the case when all block variances are equal to one, the family of polynomials of 
Corollary 5.2 gives all limit mixed moments in our random symmetric block model 
expressed in terms of asymptotic dimensions. If D is singular, then certain limit mixed 
moments of Theorem 5.1, and thus also the corresponding polynomials, may vanish. 
We close this section with giving the necessary and sufficient conditions for them to be 
non-trivial. This result shows that the situation studied in Lemma 4.2 is essentially the 
only one which is of interest in the case when at least one symmetric block in the given 
moment is not balanced. 

Corollary 5.3. Under the assumptions of Theorem 5.1, 

lim r q (n)(T pi 

if and only if the following conditions hold: 

(1) there are no evanescent blocks in the above moment, 

(2) all unbalanced blocks are of the form T p q (u, n) for some p ^ q, 
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(3) J\fCl 2 m;q {(m, Ml), ... , (tUm, Um)) # 0. 

Proof. By Theorem 5.1, it suffices to check when 

M m = y q (u Puqi ( Ul ) . ..u} Pm , qrn (u m )) # 0. 

If the three conditions hold, then M m is of the form given by Lemma 4.2 and since 
each bg(n, f) contributing to it is a positive number, we get M m ^ 0. Conversely, if 
T Pitqi (u,rii) is evanescent for some i e [m], then the corresponding u PiSi {u^ is trivial 
and thus M m = 0. Next, if T Puqi (u,rii) is unbalanced for some i e [m], where g ^ 
{pi, qi}, then, assuming (without loss of generality) that the corresponding symmetrized 
Gaussian operator is of the form uj puqi (ui) = oo Pi , qi (ui) (and thus d qi = 0) and taking the 
largest such i, we observe that it must act on Q q . or on a vector of the form e q ^ t (u) (x)w 
in order to give a non-zero result. However, by assumption q i ^ q and thus Q Qi ^ Q q 
and e qi! t(u) could only be created by u qi}t (u), which vanishes since d qi = 0. Thus, each 
unbalanced block has to be of the form T p ^ q (u, n) for some p s [r] and u e 11. Finally, 
if MCX 2 m q ((iDi, ui), . . . , (zv m , u m )) = 0, then M m = by Lemma 4.2. This completes 
the proof. ■ 

Example 5.4. If v P:q = v qtP = 1, then the moments studied in Example 5.2 reduce to 
polynomials in asymptotic dimensions d p , d q . For instance, we obtain 

%(tf, q ) = dp + H d i + d P d l 
by counting blocks of odd and even depths in the following set of partitions: 



nnn nlnl mm Inn 

The moment ty p (u pq ) is obtained from the above by interchanging d p and d q . Thus, 

= 2d p dl + 6d 2 d 2 q + 2d*d q 

since ^ = ^] . dj^frj and tyj(Q pq ) = for j $ {p, q}. Moreover, the moment under ^ q 
does not vanish even if d p > and d q = 0, which corresponds to the situation in which 
block u Ptq is unbalanced. Note that the conditions of Corollary 5.3 are satisfied. In 
particular, the first partition in the picture given above satisfies the third condition and 
contributes d p . 

In particular, if U is finite, we can consider distributions of the sums of independent 
off-diagonal and unbalanced symmetric blocks 

T p , q (n) = J]T p>q (u,n) 
ueti 

under the partial trace r q (n) associated with d q = and d p ^ 0. Under more general 
assumptions, distributions of such matrices were studied by Benaych-Georges by means 
of the so-called rectangular convolution [6]. Using Theorem 5.1, we can see that blocks 
T Pt q(u,n) behave asymptotically under T q (n) as operators uj p , q (u) under the state ty q . 
This gives random matrix models corresponding to arrays (ou Piq (u)) in which 

^p, q {u) 7^ and u qtP (u) = 

for fixed q and any p q and us U. Of course, if d q = 0, then the diagonal symmetric 
blocks T q<q (u, n) become zero as n — > go for all u e 11 under any partial trace. 



n 
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Remark 5.2. This fact allows us to give a random matrix model for non-trivial 
off-diagonal operators u> Ptq (u) in the case of arbitrary dimension matrix. Namely, if 
Tp, q (u, n) consists of a block S Ptq (u,n) and of its adjoint, where S Ptq (u,n) is expressed 
in terms of its columns as 

S Piq (u,n) = (Kxfari) ... K riq (u,n)), 

then uj Pjq (u) is the limit under r q (n) as n — » go of the symmetric subblock of T p ^ q (u,n) 
which consists of the column Kj(u,n) and of its adjoint for arbitrary j. In the case 
when all block variances are normalized by n, it does not mean that one can construct 
a random matrix model for any matrix (u) p>q (u)) since if one entry of this matrix is 
zero, all entries in the same row must be zero. However, if we consider more general 
normalizations, this goal can probably be achieved. 

6. Double Convolution Formula 

The limit distributions of sums of independent Gaussian square random matrices 
with complex identically distributed variables are given by the free additive convolution 
[32,33]. In the case when the variables are block-identically distributed and 11 is finite, 
in order to express the limit distributions of the sums 

Y(n) = Y i Y(u,n) 

ueu 

under partial traces in a convolution form, we need to generalize Voiculescu's result and 
introduce the convolution associated with matricial freeness. 

This convolution, called the matricially free convolution, reminds the strongly matri- 
cially free convolution 

EB P , q H P , q = distribution of ^a Pi9 

pa 

in the state tp, where (a Pjq ) is an r x r array of strongly matricially free random variables 
with respect to an array of states built from a distinguished state ip and its conjugate 
states, and (/x P)9 ) is the corresponding array of distributions. 

The concept of strong matricial freeness and the associated convolution was studied 
in detail in [22,24]. We have shown that it unifies the main types of noncommutative 
independence in the sense that various noncommutative random variables can be de- 
composed in terms of strongly matricially free ones. One of the consequences of this 
unification is that both free and boolean convolutions can be expressed in terms of the 
strongly matricially free convolution of an array of distributions. 

In particular, if the array (a Pjq ) is square and of dimension at least two and the 
corresponding distributions (fJ, Ptq ) are row-identical, then 

EBp.g/ip.q =|llffl/l2E0...EB/ir, 

where \x PA = \x p for any p, q, which justifies the notation since in this case the convolution 
EB p ,q can be viewed as a decomposition of the free convolution. However, if the array 
(a Piq ) is diagonal, then 

where fi p>q = \x p whenever p = q. It should be remarked that the strongly matricially 
free convolution is not just the free convolution taken over an array of indices. These 



LIMIT DISTRIBUTIONS OF RANDOM MATRICES 35 

properties of the strongly matricially free convolution will be used in this section. In 
the meantime, let us introduce another convolution of matricial type. 

Definition 6.1. Let (a P:q ) be an array of variables which are matricially free with 
respect to the array (<p p , q ), where ip p>q = ip q for any (p, q) e J, with distributions (/i p , g ), 
respectively. By the matricially free convolution of (// p , g ) we shall understand 



Sp,g/ip, g = distribution of ^a Pi 



i 

p,q 



in the state ip = V d q <p q , where d q ^ for any q, where we supress its dependence on 
ip or on the dimension matrix D. 

The most important examples of matricially free convolutions are the distributions 
of 'independent' Gaussian pseudomatrices [23] 

u(u) = Yj^pM 
p,q 

which generalize free Gaussian operators, and the distributions of their sums 

u = Y u ( u ) 

ueti 

in the state \&, where the associated arrays of distributions are matricial semicircle laws. 

Example 6.1. By [23, Proposition 4.2], where we have just one matrix of matricially 
free Gaussian operators, the distribution of oo = 2 P q ^p,q in the state \1> = Yj q d q ^ q is 
given by /i = \*} p , q n P , q , where 

_ J (T q , q iip = q 
»™ \ n p , q iip^q ' 

which coresponds to the simplest case of Proposition 2.3, when the set 11 consists of 
one element. 

In order to express the distributions of the sum of matricially free Gaussian operators 
in a convolution form in the general case, we shall introduce a new notation. If (// P)9 ) 
and (v p ,q) are arrays of distributions indexed by J, we shall denote 



m f Vq,q ffl Vq,q Xp = q 



where /zEBz/ and fi w v denote the free convolution and the boolean convolution, respec- 
tively. It can be seen that this convolution is pointwise associative, i.e. associative for 
fixed (p, q) e J7 since both free and boolean convolutions are associative. 

Lemma 6.1. Under the assumptions of Proposition 2.3, the ^ -distribution of the sum 
of collective Gaussian operators oo = ^ p oj va is given by 

V = H P , g (£t Pi? (l) EBj • • • EB? H P ,q{m)) 

where the arrays (fj, p>q (u)), us U = [to], are given by Corollary 2.1. 

Proof. We will first prove that the family 
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is free (boolean independent) with respect to ip p>q = for any fixed diagonal (off- 
diagonal) (p, q) e J . Let M Pt q{v) be the subspace of N PA spanned by tensors which 
begin with e p>q (u) and let s Pjq (u) be the associated orthogonal projection. First, we 
will treat the off-diagonal case. Denote, for simplicity, p{u) = p Piq (u), p*{u) = p pq (u), 
s(u) = s Ptq (u) and l(u) = l Ptq (u), where l p , q (u) is the orthogonal projection onto 

cn,eAf M (ii)e0Af 5it , 

k 

thus note that l p , q (u) is slightly smaller than of Definition 2.5, and let 

A{u) := alg(p(u), p*(u), l(u)) 

for any u e U. We will show that 

■V q ( ai ...a n ) = ^ q (a 1 )...^f q (a n ) 

for any e A(uk), where u±, . . . , u n e U. First, let us check how a n acts onto fl q . Since 
a n is a linear combination of p(u n ), p*(u n ), s(u n ) and l(u n ), we clearly have 

a n Q q = ^q(a n )Q q + ^e M {u n ) 

for some 7 e C. However, ill ^ k, each element of A(u{) kills all vectors which begin 
with e Pt g(uk) by Proposition 2.1. Therefore, 

<ai . . . a n fi 9 , ttq) = *,(ai) . . . * g (a n ) 

which shows boolean independence of subalgebras .4.(1), . . . ,A(m) with respect to 
Therefore, the associated Gaussian operators are also boolean independent. Now, let 
us consider the diagonal case and fix (q,q) e J. Again, denote p(u) = pq <q (u) and 
p*(u) = pq^ q {u) for any u. In this case, we take the same unit l q>q for each u, and we 
denote 

A{u) = alg(p(u), p*{u)) 

where we remark that we do not need to put the unit \ qA among the generators since 
Pq,q( u )pq,q( u ) = bq t q(u)lq jq and we assume that b qq (u) > 0. Using the relations of 
Proposition 2.2, we can write every element of A(u) n Ker\l/ g in the form of a linear 
combination of 

p p (u)p* q (u), where p + q > 0. 

However, since p*{u) kills Vt q and, moreover, p*{u)p{t) = for any u ^ t, in order to 
prove the freeness condition 

vf^ai ...a n ) = 
for any e A(iik) n Ker^ q , it suffices to observe that 

Vq(p P ^U 1 )...p P ™(u m ))=0. 

Finally, matricial freeness of the array of algebras A PtQ = alg(p Pjq , p p q , \ PA ) with respect 
to ((p Pl q) follows from Proposition 2.3, and that entails matricial freeness of the array 
(uJ p ,q). This completes the proof. ■ 

Definition 6.2. Let (ajj) be a symmetric array of variables which are symmetrically 
matricially free with respect to the array (ipij), where tpij = ipj for any i,j, with 
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distributions (fiij), respectively. By the symmetrically matricially free convolution of 
(fJ>ij) we will understand 

0ijA*i,j = distribution of ^ ciij 

in the state ip = ^] ■ djtpj, where we supress its dependence on the dimension matrix D. 

Lemma 6.2. Under the assumptions of Proposition 3.5, the ^/-distribution of the sum 
Q = XrtT=i £( u ) i s given by 

li = P ,,(ff M (l) ffl • • • ffl ff M H) 
where (a Ptq (u)) is the array of semicircle laws of Corollary 3.1 for any ue U. 

Proof. By Proposition 3.5, the array of variables 

m 
u=l 

is symmetrically matricially free with respect to the array (^ Ptq ) defined by . . . , \I/ r . 
Moreover, the distribution of Q Piq (u) in the state ty Ptq is the semicircle law of radius 
2a PA {u) for any (p, q) e J7 and u e 11. Therefore, it suffices to show that for any given 
(j), q) e J, the family of variables 

{w M (l),..,w M (m)( 

is free with respect to ^/ q . The proof for the case when p = q is the same as in Lemma 
6.1 since Q qA {u) = ui q ^ q {u) for any u. Let us fix (p,q) e J7 and denote p(u) = p p , q (u) 
and p*(u) = p£ q (u) for any u. Let 

A(u) = alg(p(u),p*(u)) 

where, as in the diagonal case, we do not need to include the unit l PtQ among the 
generators since pt Ju)p P:q (u) = b Ptq (u)l p , q and b p ^ q (u) > 0. Using the commutation 

relations of Proposition 3.4, we can write every element of A{u) n Ker\l/ g in the form 
of a linear combination of 

P j (m)p*' j (-u), where i + j > 0. 

However, since p*{u) kills Q q and, moreover, p*{u)p{t) = for any u t, in order to 
prove the freeness condition 

■V q (a 1 ...a n ) = 
for any e A(uk) n Ker^ q , it suffices to observe that 

y q (p p ^u 1 )...^(u m )) = 0, 
which completes the proof. ■ 

We choose Lemma 6.1 to find a convolution formula for the limit law of the sum of 
independent Hermitian Gaussian random matrices. In the identically distributed case, 
we obtain the free convolution of semicircle distributions [32]. In our generalization of 
that result, the free convolution is replaced by the strongly matricially free convolution. 
In addition, we obtain the matricially free convolution which tells us how to convolve 
asymptotic distributions of different blocks. These distributions are computed with 
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respect to partial traces or with respect to their convex linear combination given by the 
trace. We can call this formula the Double Convolution Formula since it involves both 
the matricially free convolution and the strongly matricially free one. 

Equivalently, one can work with the symmetrically matricially free convolution. In 
particular, one can use Lemma 6.2 to prove a similar convolution formula to that of 
Theorem 6.1. However, an analytic approach to the symmetrically matricially free 
convolution is more complicated and will not be treated in this paper. For that reason, 
we restrict ourselves to the case of the matricially free convolution in the theorem given 
below as well as in Section 8, where we discuss some analytic results. 

We shall use the square roots of all entries of B(u), namely 

a P,g( u ) = 

for any (p, q) e J and u e 11, with the corresponding matrices denoted A{u) = (a PA (u)). 
Moreover, in order to express the result in terms of the strongly matricially free convo- 
lution, we shall use pairs (r, s) e 11 x 11 instead of u e 11, and for that purpose we set 
a g ,q(r, s) = a q>g (r) and a p>q (r, r) = a p>q (r) for any p # q and any r,se 11. 

Theorem 6.1. Under the assumptions of Theorem 5.1, the limit law of the sum of 

independent Hermitian Gaussian random matrices 

Y(n) = Yj Y (u,n) 

ueu 

under the trace takes the convolution form 

n = a,- 8^(0 

where j e 3 and £ e Zl x 11 and 

(1) [J,j(£) = o}(0 for any fi e 11 x 11 whenever j £ J is diagonal 

(2) = for any diagonal £e 11 x 11 whenever j e J7 is off-diagonal, 

where <7,-(£) is the semicircle law of radius 2a 7 -(£) and is the Bernoulli law con- 

centrated at +aij(k). 

Proof. We have 

Y m {n)= Yj Yj T Puqi( u un) ...T Pmt q m (u m ,ri) 

ui,...,u m e'Upi,qi,...,Pm,qm. 

and thus, using Theorem 5.1, we obtain 

lim r(n)(y"») = ^(u m ) 

for any m s N, where 

p,q p,q ueu 

Now, by Lemma 6.1, the ^/-distribution of uj is given by 

\*\p,q{^p,q{l) Si ... Si Vp,q(m)), 

where (n P , q (u)) is the array of semicircle and Bernoulli laws of Corollary 2.1 for any 
u e 11 and any p, q. It remains to decompose the free and boolean convolutions in 
terms of the strongly matricially free convolutions [24, Theorem 3.1]. Using the above 
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notation, where the pair (p, q) is already reserved for the matricially free convolution, 
we can write these decompositions in the form 



where 'J q , q (r, s) = H q , q (r) for any 1 s; r, s < m and any (q, q) e J and 7 P , g (r, r) = fJ, p>g (r) 
for any 1 < r < m and J p , q (r, s) = 5q when r ^ s, for any off-diagonal (p, q) e J/. This 



which corresponds to the matricially free convolution of semicircle and Bernoulli 
laws which form one matricial semicircle law. 

Proof. The first statement holds since in the case when r = 1, all matrices have 
one diagonal symmetric block whose complex entries are identically distributed and for 
r = 1 the strongly matricially free convolution over £ reduces to the free convolution. 
The second statement restates the result of [23] on the distribution of one Gaussian 
random pseudomatrix oj = Y\ nn u) vg in the state \& '. ■ 

In the general case, when both convolutions are present, we can interpret their origins 
as follows: 

(a) convolution EB^ describes the asymptotic relation between symmetric blocks hav- 
ing the same position in independent matrices, 

(b) convolution [*L describes the asymptotic relation between symmetric blocks hav- 
ing different positions in a given matrix, 

where by an independent matrix we understand an independent Hermitian Gaussian 
random matrix. 

Remark 6.1. It should be noted that the limit distribution of Y(n) under the trace r(n) 
given by the Double Convolution Formula depends on the dimension matrix D. In fact, 
D uniquely determines the convolution [*]. However, it is not hard to see that a similar 
formula holds also for limit distributions under the partial traces r q (n), where q e [r]. In 
fact, if d q = 0, the moments under ty q do not contribute to the asymptotic distributions 
under the trace whereas they may be non-trivial as Example 5.4 demonstrates (see also 
Section 8). To avoid ambiguity, we prefer to speak of the convolution under the state 
\& or under the state ^/ q , where qs [r], instead of describing \*\ in terms of D. 

Remark 6.2. Both boolean and free convolutions are associative, and thus we can 
define an associative convolution of arrays [/x] = {fJ- P , q ), \y\ = (v p , q ), where (p,q) s J7, of 
probability measures on the real line 



n q , q {i)m...m^ q , q {m) 

Hp, q {l) H P , q (m) 



Sr,slq, q (r, S) 
Rr,s7fp,q(r, S) 



completes the proof. 




M s M = [ 7 ] 
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where 

7p,g = EB? v v,q 

for any (p, q) e J7, where we use the notation El as a generalization of the free convolution 
to the family of arrays. Using this convolution, one can rephrase the Double Convolution 
Formula as follows. The sum 

m 

Y(n) = J]y(n,w) 

of independent Hermitian Gaussian random matrices converges in moments to the pseu- 
domatrix u, meaning that the moments of Y(n) under r q {n) converge to the moments 
of uj under ty q for any q e [r], and thus converge to the corresponding moments of the 
free convolution of matricial semicircle laws 

[a] = [a(l)]EH2)]EI...El[a(m)]. 

This gives convergence in moments for any convex linear combination of the cosnidered 
states. One can also define the associated notion of independence for arrays of noncom- 
mutative probability spaces, but we will not discuss it here. 

7. Asymptotic monotone independence and s-freeness 

The Double Convolution Formula seems to suggest that other convolutions decom- 
posable in terms of the strongly matricially free convolution might be obtained as limit 
distributions of random matrix models using the framework of matricial freeness. In 
the case of block-identically distributed arrays, we were able to obtain only free and 
boolean convolutions as the limit distributions of the diagonal and off-diagonal entries, 
respectively, of the sum of Hermitian Gaussian random matrices. 

In this section we give two other examples, the monotone and s-free convolutions. 
However, in order to do that, one has to recall that they arise in the context of the 
decompositions of the free convolution in terms of other convolutions [21] and that they 
can be obtained from the strongly matricially free convolution [22]. Therefore, using 
these decompositions, it is natural to decompose independent Hermitian Gaussian ma- 
trices with identically distributed entries in terms of symmetric blocks and then suitably 
identify the right pieces of these decompositions with random variables independent in 
the required sense. 

We will restrict our attention to the case of two independent Hermitian Gaussian 
random matrices decomposed in terms of their three symmetric blocks as 

Y(u,n) = T M'K n )> 

where u e {1,2}. Choosing appropriate symmetric blocks of these matrices, we will 
construct pairs of independent Hermitian Gaussian random matrices {X(l, n), X(2, n)} 
which are asymptotically boolean independent, monotone independent or s-free under 
one of the partial traces. Thus, the distributions of the sums 

X{n) = X(l,n)+X(2,n) 

will be asymptotically boolean, monotone or s-free convolutions of the asymptotic dis- 
tributions of the summands. The boolean case follows immediately from Lemma 6.1, 
but the other two are slightly more involved. 
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In view of Theorem 5.1, if we are able to construct families of random variables which 
are independent in some sense from a family of arrays of matricially free Gaussian 
operators, a corresponding random matrix model is close at hand. The definitions of 
boolean independence, monotone independence and s-freeness can be found in [22]. Let 
us only remark that in the last two cases the order in which the variables appear is 
relevant. 

Proposition 7.1. If U = {1,2} , r = 2 and b P:q (u) = 1 for any p,q,us [2], then 

(1) the pair {a>i i2 (l), ^1^(2)} is boolean independent with respect to ^2, 

(2) the pair {^^(l), oo\^{2) + u\^{2)} is monotone independent with respect to \I/ 2 , 

(3) the pair {wi i2 (l) + wi,i(l), wi,i(2)} is s-free with respect to (^^i)- 

Proof. (1) follows from Lemma 6.1. To prove (2), denote 

a = wi j2 (l), b = wi j2 (2) + a;i,i(2). 

We need to show that the pair {a, b} (in that order) is monotone independent w.r.t. 
* 2 , i-e. 

^2(^1016102^2) = ^ 2 (bi)^2(wia 1 a 2 W2) 

for any 01,02 e C[a, li], b\ e C[6, 12], where li = li )2 and I2 = 1 and wi,w 2 are 
arbitrary elements of C{a,b, li, 1 2 ). It suffices to consider the action of a and b onto 
their invariant subspace in Ai of the form 

m' = cn 2 © 0F(2) ® n(i)) © (^(2) ® u{2)) 

where F{2) = J r (Cei i i(2)) with the vacuum vector Q and H{u) = Cei t2 (u) for u e {1, 2}, 
where we identify Vt (x) ei )2 (ii) with ei ;2 (u). Now, the range of any polynomial in a is 
Cn 2 © W) since 

fe n f ^2 if & even 
tt 2 = ( ei )2 (l) if k odd 

and li^2 = ^2 ; 1161,2(1) = ei ;2 (l). Therefore, it suffices to compute the action of any 
polynomial in b onto Q 2 and ei j2 (l). Now, the action of powers of b onto Q 2 and onto 
ei j2 (l) is the same as the action of the free Gaussian operator onto the vacuum vector 
in the free Fock space. Namely, we have 

b 2k Q 2 = C k Q 2 mod (M! Q CQ 2 ) 

and 

b 2k e 1>2 (l) = C fc ei, 2 (l) mod (M' © (Cfi 2 © U{1)) 
where Ck is the kth. Catalan number and 

b 2k - 1 Q 2 = b 2k - l e lt2 {l) = mod (M' (CQ 2 U{\)) 

for any k s N. Thus ^ 2 (6 2fc ) = C fc and, moreover, since M' (Cf2 2 "%(!)) c Kera, 
the required condition for monotone independence holds if 61 is a positive power of b. 
It is easy to see that it also holds if b± = 1 2 , which completes the proof of (2). To prove 
(3), consider 

a = wi, 2 (l) +wi > i(l), b = u 1A (2) 
and their invariant subspace in Ai of the form 

M" = cn 2 0(^(1 ©2) ®H(1)) 
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where J(10 2) = 7"(Cei,i(l) Ce lj i(2)). Observe that M" is isomorphic to 

7"i, 2 = cn e 7° (7j 7°) (7? J=l 7°) . . . 

where 7*° = J r (Ce- ) ) C£j and the isomorphism r is given by r(f2 2 ) = ^ and 

T(ei,i(ji) ... ei,i(ifc_i) ei )2 (l)) = e ix . . . e ifc _ 1 ei 

for any ji, . . . , jk-i e {1, 2} and fceN. The space J 7 !^ is the s-free product of two free 
Fock spaces, and (J-2,^2), the subspace of their free product (J-i,£i) * (J-2,^2) 

(for the definition of the s-free product of Hilbert spaces and of the s-free convolution 
describing the free subordination property, see [21]). Now, observe that the action of 
a onto M. can be identified with the action of the canonical free Gaussian u(ei) re- 
stricted to 71,2 • Similarly, the action of b onto M." can be identified with the action of 
cu(e 2 ) restricted to T\, 2 . This proves that the pair {a,b} (in that order) is s-free with 
respect to the pair of states (^ 2 , ^1), which gives (3). ■ 

We are ready to state certain results concerning asymptotic properties of independent 
Gaussian random matrices X(l,n), X(2,n) built from at most three symmetric blocks 
of Y(l, n), Y(2, n), respectively. We assume that one asymptotic dimension is equal to 
zero, thus one block is balanced, one is unbalanced and one is evanescent. 

Theorem 7.1. Under the assumptions of Theorem 5.1, let U = {1,2}, r = 2, d 2 = 
and d 1 = 1 and let all block variances to be equal to one. Then 

(1) the pair {Ti j2 (l, n), 7i j2 (2, n)} is asymptotically boolean independent with respect 
to T 2 {n), 

(2) the pair {Ti j2 (l, n), Y(2, n)} is asymptotically monotone independent with respect 
to T 2 (n), 

(3) the pair {Y(l, n), 71^(2, n)} is asymptotically s-free with respect to (r 2 (n), Ti(n)). 

Proof. Since d\ = 1 and d 2 = 0, blocks T^\{u, n), Ti t2 (u, n) and T 2t2 (u,n) are bal- 
anced, unbalanced and evanescent, respectively. Therefore, by Theorem 5.1, we have 
convergence as n — > 00 

Ti, 2 (u,n) -> u 1)2 {u) 
T 1A (u,n) -> u 1A (u) 
T 2 , 2 (u,n) -> 

in the sense of moments under r 2 (n) for u e {1,2}, where the limit moments are com- 
puted in the state ^ 2 . By Proposition 7.1, the proof is completed. ■ 



8. Catalan matrices 

In general, finding explicit analytic formulas for the limit distributions given by the 
Double Convolution Formula is difficult, if at all possible. However, one can prove 
existence of a matrix-valued reciprocal Cauchy transform whose coefficients are matri- 
cial analogs of the Catalan numbers. Moreover, certain formulas for the scalar-valued 
Cauchy transforms of these distributions or their reciprocals can also be dervied. 
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Let us recall that the convolutions of Theorem 6.1 are distributions of the operator 
of the form 



where each uj va is a collective Gaussian operator defined as in the proof of Theorem 
6.1. In view of Lemma 6.1, it is easy to determine the distributions of operators u P)Q in 
the states ty pq = respectively, and thus the distribution in the state ty. 

Proposition 8.1. Each diagonal collective Gaussian operator u qA has the semicircle 
distribution of radius 



in the state ^/ q for any q e [r] . Each off-diagonal collective Gaussian operator uj p ^ q has 
the Bernoulli distribution concentrated at 



in the state ty q for any p q. 

Proof. Each diagonal operator uj qA {u) has the semicircle distribution of radius 2a q , q {u) 
for any q e [r] and u e U by Proposition 2.3. The family {u qjq (u),u e U) is free w.r.t. 
ty q for any q by Lemma 6.1. Therefore, the corresponding sum u qA has the semicircle 
distribution of radius 1ct qA given above. In turn, since each off-diagonal operator u PA (u) 
has the Bernoulli distribution concentrated at +a Ptq (u) and the family {ui Ptq (u),u s 11} 
is boolean independent w.r.t. ^/ q for any p ^ q, the corresponding sum oj pa has the 
Bernoulli distribution concentrated at +a Ptq given above. ■ 

Therefore, the array (oo p , q ) of collective Gaussian operators is of the same type as in 
the case of one array. This implies that the results on the Cauchy transforms of the limit 
distributions of random pseudomatrices (sums of matricially free random variables) 
studied in [22] can be used in our random matrix context. In fact, we showed in [23] 
that symmetric blocks of a Hermitian Gaussian random matrix of Theorem 5.1 have 
the same limit joint distributions as blocks of a random pseudomatrix. Therefore, 
the distributions given by the continued fractions of Theorem 8.1 coincide with the 
matricially free convolutions of the distributions of (oo Ptq ). 

Theorem 8.1. The Cauchy transform of the distribution of the operator u in the state 
^/ q , where ge [r], takes the form 




in the state 



r 





i 
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where b m = Yj u eu a p,q( u ) f or an U P,q e [ r ] an d where the continued fraction converges 
uniformly on the compact subsets ofC + to a probability measure Li q on the real line with 
compact support. 

Proof. We have shown in [23, Theorem 7.1] that mixed moments of blocks of pseudo- 
matrices, converge under partial traces ip q (n) to the mixed moments of the matricially 
free Gaussian operators under ty q , respectively. These operators have semicircle or 
Bernoulli distributions and are matricially free by [23, Proposition 4.2], or by the gen- 
eralized version of this result given in Proposition 2.3 of this paper. Therefore, since the 
collective Gaussian operators u Ptq have semicircle or Bernoulli distributions by Propo- 
sition 8.1 and the array (oj p>q ) is matricially free with respect to (&p, q ) by Proposition 
2.3, we can treat them as the limit operatorial realizations of blocks of certain random 
pseudomatricies under partial traces. We showed in [22, Theorem 7.1] that the distri- 
butions of such blocks under partial traces converge weakly to the measures /xi, . . . , \x r 
on the real line whose Cauchy transforms are of the form given above. This proves our 
assertion. ■ 



Corollary 8.1. The measures of Theorem 8.1 take the boolean convolution form 

l*q = fil,q ^ fl2, q tb . . . tb 

for any q e [r], where (fJ>ij) is the array of probability measures on the real line whose 
K-transforms are defined by continued fractions 

K P,q( Z ) = 7 

v b k,p 

z - s fe — 



bj,k 



3 z- 



which converge on the compact subsets of C + , where K^(z) = z — F^(z) and F^(z) is 
the reciprocal Cauchy transform of fx. 

Proof. This is an immediate consequence of Theorem 8.1 and [22, Lemma 7.1]. ■ 

Finding an analytical formula for G^ q in the general case does not seem possible. 
Let us present two situations discussed in [22] in which certain interesting convolution 
formulas can be derived. From these formulas one can easily obtain the transforms. 

Corollary 8.2. If A = (a Pjq ) is a square r-dimensional matrix with identical positive 

entries a PA = ^Yj u &u a p,q( u ) = a p ^ n P'^ 1 row ' then 

^ q = o~ ai H a a2 EE . . . EE a ar 

for each q e [r], where cr ap is the semicircle distribution of radius 2a p for p e [r] . 

Proof. This is an immediate consequence of Theorem 8.1 and [22, Proposition 8.2]. ■ 

Example 8.1. Using Corollary 8.2, we can give a simple example of the matricially 
free convolution. Let us consider the case of one array with all variances equal to one. 
We convolve the array (/i p , g ), where fx qjq is the semicircle law of radius 2-\fd~ q and /j, P:q is 
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the Bernoulli law concentrated at + *Jd~ v f° r an y P ^ Q- The matricially free convolution 
of Theorem 6.1 takes the form 

r 
9=1 

since fx q = o~\ for any q e [r] by Corollary 8.2. In fact, the free convolution of semicircle 
laws with radii 2\fd\, . . . , 2\fd~ r is the semicricle law with radius one, where we use 
d\ + . . . + d r = 1. A similar result can be obtained if we convolve a finite family of 
arrays. 

In the case of two-dimensional square arrays, we can express the distributions fii, fx<i 
in terms of s-free additive convolutions. These convolutions were introduced in [21] in 
the context of the subordination property of the free additive convolution of proba- 
bility measures on the real line .Mr discovered by Voiculescu [34] and generalized by 
Biane [8]. For details on the s-free convolution and the convolution formulation of the 
subordination property given below, see [22]. 

Theorem 8.2. For any fx, v e .Mr there exists a unique v □ fx e .Mr such that 
called the s-free convolution of fx and v. Moreover, it holds that 
where fx i> v is the monotone convolution of fx and v. 

Remark 8.1. This is a convolution formulation of the subordination property of the 
free additive convolution, usually written in terms of the reciprocal Cauchy transforms 
of probability measures on the real line as 

where Fi(z) and F 2 (z) are the reciprocal Cauchy transforms of the unique probabil- 
ity measures on the real line which satisfy the above equations. In our convolution 
formulation, these measures are denoted by v □ fx and fx[B^, respectively. 

In order to write convolution formulas for the distributions of u under ^ and ^2, we 
shall use boolean compressions of semicircle laws. For any probability measure fx on the 
real line, we denote by T t fx its boolean compression defined in terms of K-transforms as 

K Ttll = tKfj,, 

where t ^ 0. We allow t = 0, in which case Tqlx = Sq. We shall use two-parameter 
boolean compressions of semicircle distributions, = T t a a for t = (f3/a) 2 , where 

(2a 2 - (3 2 )z - (3 2 ^z 2 - 4a 2 
^ Z >~ (2a 2 - 2/3 2 )z 2 + 2/3 4 

are their Cauchy transforms, where the branch of the square root is chosen so that 
\/z 2 - 4a 2 > if z e R and z e (2a, 00). 

In view of the decomposition of the free convolution in terms of s-free convolutions 
of Theorem 8.2, the result given below shows that the distributions of u> under ^ and 
^2 are deformations of the free convolution of semicircle laws. 



46 



R. LENCZEWSKI 



Proposition 8.2. If (uj Pjq ) is a square two-dimensional array with no zero operators, 
then the distributions of u under ^ and \& 2 ; respectively, take the form 

^2 = ri /a (ft >7 □ ?a i/3 ) w (? a>/3 □ ft >7 ) 

respectively, where a = bi tl , f3 = b lj2 , 7 = &2,i; <5 = ^2,2; an d t = (j3/a) 2 , s = ('"f/S) 2 . 

Proof. These formulas were proven in [22, Theorem 8.1] on the basis of the continued 
fraction representation of the Cauchy transforms of Theorem 8.1. ■ 

If some of the matricially free Gaussian operators in the two-dimensional array vanish, 
which can be taken into account by setting the corresponding distributions to be S , 
their matricially free convolutions can also be retrieved from [22]. Let us analyze one 
special case in more detail. 

Example 8.2. The limit distribution of Example 5.2, derived from Theorem 5.1 in a 
purely combinatorial fashion, can be obtained using convolutions. In this case, \i is the 
matricially free convolution of two Bernoulli distributions, fi Ptq = k pa (concentrated at 
+a p , q ) and fj, qp = K qp (concentrated at +a qtP ) and thus, by [22, Theorem 7.1], it is 
equal to their s-free convolution 

which is the distribution $ corresponding to a = b VA and b = b qjP . Moreover, it can be 
seen from Theorem 6.1 that a similar result holds for the sum 

T P,q( n ) = Yj T P^ U ' n ^ 

ueu 

since we first compute the boolean convolution of Bernoulli laws to get 

K P,q = l+J K P,q( U )' 

ueu 

the Bernoulli law concentrated at +a Pi9 of Proposition 8.1. In a similar way we obtain 
K q>p . Then, using matricial freeness, we obtain again the s-free convolution of two collec- 
tive Bernoulli laws, K PtQ and K QtP . In other words, it suffices to replace the pair (b Piq , b qtP ) 
by the pair (XLeu W,q( u )y H U eu ^q,p( u ))- For other matricially free convolutions, see [22, 
Corollary 8.2]. 

Finally, let us present a matricial approach to the limit distributions of Theorem 8.1. 
This approach leads to matricial analogs of Catalan numbers which are matrix-valued 
coefficients in the Laurent series for the matrix consisting of diagonal K-transforms. 

Note that the Cauchy transforms of Theorem 8.1 can be expressed in terms of the 
K-transforms K Pjq (z) of Corollary 8.1. Let us suppose that b qA > for any q e [r]. 
In that case, these transforms can be reduced to the diagonal ones. For that reason, 
introduce diagonal matrices 

C(z) = dmg(K 1A (z),...,K rtr (z)) 
B = diag(6i ) i,...,6 r>r ) 

as well as the matrix 

T = (t Ptq ) e M r (C), where t p>q = b p Jb p , p 
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for any p,q s [r] . Moreover, let 

<D : Af r (C) — M r (C) 

be the mapping given by 



2) 



(A) = diag ^a;,!, . . . , J] <V 



vi=l i=l 



for any A = (a^) e M r (C). 

Using these notations, we can find a formula for C(z), from which the K-transforms 
of n q and of n = Yj q dql^q can be retrieved. The coefficients of the matrix- valued series 
obtained in this fashion play the role of matricial analogs of Catalan numbers called 
Catalan matrices. 

Theorem 8.3. If b qq > for any q e [r], then C(z) assumes the form of a formal 
Laurent series 

oo 

C(z) = J] C n z- 2n - 1 

?1 = 

where (C n ) is a sequence of diagonal matrices satisfying the recurrence formula 

C n = J] Di&TCj) 

i+j=n—l 

for any natural n, with C = B. The series C(z) converges in the supremum norm for 
\z\ > 2^/r || T (HI B || . 



Proof. By Corollary 8.1, we have 



Kp,q\Z) tp^qKp^p^Z^J 



since K Pjq (z) has the same denominator as K PtP (z) for any p, q (of course, t g>q = 1 for 
any q e [r]). Moreover, 



K q , q (z) = 



J q,q 



z 



2p=l tp,qKp,p( Z ) 

for any qs [r] . This leads to the equation 

<D{C{z)TC{z)) - zC{z) + B = 

Since each K q q (z) is the K-transform of some probability measure on the real line by 
Corollary 8.1, the matrix C(z) is of the desired form 



z 2n - 1 



n=0 



where C n is a constant diagonal matrix for each n. Substituting this series into the above 
equation and computing its coefficients, we obtain the desired recurrence formula. As 
for the convergence of the series C(z), note that the recurrence formula is of the same 
type as for Catalan numbers and thus the expression for C n has C n terms. Using this 
recurrence, we can estimate the supremum norm of each C n : 



C n UC n r n II T HI 5 



in+l 
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and, denoting the right hand side by a n , we obtain 

a n 1 



lim , ... 

ti^od a n+1 4r II T \\ \\ B 



and thus, using D'Alembert's test for convergence, we obtain convergence for 



\z\ > 2^/r || T HI! B ||, 



which completes the proof. 



Remark 8.2. In order to see the connection between Catalan matrices and Catalan 
numbers in a more explicit way, let us introduce the bilinear mapping 

<.,.>: M r (C) x M r (C) -> M r (C) 

by the formula 

(A 1 ,A 2 )=>D(A l TA 2 ). 
Then, we can write the first few matrices C n in the form 
Ci = (B,B) 

C 2 = (B,(B,B)) + ((B,B),B) 

C 3 = (B,(B,(B,B))) + (B,((B,B),B)) + ((B,B),(B,B)) 
+ ((B,(B,B)),B) + (((B,B),B),B), 

etc. If r = 1 and = 1, then C n is equal to the nth Catalan number C n since in 
one of the many counting problems involving Catalan numbers, C n is the number of 
different ways n + 1 factors can be completely parenthesized. In turn, if b VA = d p for 
any p, q e [r], then 

C n = diag(diC n , . . . , d r C n ) 

for any natural n, which corresponds in the random matrix context to the situation in 
which all block variances are equal to one. 



9. Non-Hermitian Gaussian random matrices 

We would like to generalize the results of Section 5 to the ensemble of independent 
non-Hermitian Gaussian random matrices. This ensemble (or, the family of its entries) 
is sometimes called the Ginibre Ensemble (the family of symmetric blocks of the Ginibre 
Ensemble can be called the Ginibre Symmetric Block Ensemble). We keep the same 
settings for blocks as in Section 5. 

We will show that the Ginibre Symmetric Block Ensemble converges in *-moments 
as n — > oo to the ensemble of non-self adjoint operators 

where 

rj P , q (u) = pp, q (2u - 1) + p* Pjq (2u) 

for (p, q) e J7, u e U = [t] and (p p , q (u)), where u e [2t], are arrays of symmetrized 
creation operators (the set J7 is assumed to be symmetric). In fact, it suffices to consider 
the case when J7 = [r] x [r] and restrict to a symmetric subset if needed. 
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In order to define t arrays (r} p , q {u)), where us [t], we need 2t arrays (p P;q (u)), where 
u s [2t\. Moreover, we shall assume that 

b P , q (2u - 1) = b PA {2u) = d p v m (u) 

for any fixed p, q, u. Our definition parallels that in the free case, where operators 

rj(u) = £(2u- 1) + t(2u) 

for us [t], which are unitarily equivalent to circular operators, are introduced [32, 
Theorem 3.3]. Here, {£(1), • • • ,£(2t)} is a family of free creation operators. 

We can now generalize the result of Theorem 5.1 to the Ginibre Symmetric Block 
Ensemble. The important difference is that we do not assume that the matrices are 
Hermitian and thus each matrix contains 2n 2 independent Gaussian random variables. 
However, in order to reduce this model to the Hermitian case, we need to assume that 
the variance matrices are symmetric. The theorem given below is a block refinement of 
that proved by Voiculescu for the Ginibre Ensemble [32, Theorem 3.3] except that we 
assume the Gaussian variables to be independent block-identically distributed (i.b.d.) 
instead of i.i.d. 

Theorem 9.1. Let (Y(u,n)) be the family of complex random matrices, where neN, 
us U = [t], and let 

Tp, q (u, n) = Y i,j( u ' n ) ® e M'( n )> 

where N PiQ := (N p x N q ) u (N g x N p ) for p, q s [r] and reN, and where 

{ReYij(u, n), ImY^(-u, n) : i,j s [n),us U] 
is a family of independent Gaussian random variables with zero mean and variances 

E((ReY hj (u,n)r) = E((lraY t ,(u,n)f) = ?eM 

for any s N p>q and us U, where the matrices (v Ptq (u)) are symmetric. Then 

limT q (n)(T p l qi ( Ul ,n) . . . T^Jum, n)) = . ■ ■ V P Z, q J u m)) 

for any e 1: ...,e m e {1,*}, u u ...,u m s U and q,p u q u . . . ,p m , q m s [r]. 

Proof The proof is similar to that in [32, Theorem 3.3]. Write each matrix in the 
form 

Y(u, n) = 2- 1 ' 2 (x(u, n) + iX(u, n)) , 

where the matrices 

X(u,n) = 2- 1/2 (Y(u,n) + Y*(u,ri)) 

X(u,n) = i2- 1/2 (Y*(u,n)-Y(u,n)) 

are Hermitian for any neN and us U. The symmetric blocks of these matrices, denoted 
(U Pi g(u, n)) and (U p>q (u,n)), respectively, will give the asymptotics of symmetric blocks 
(Tp >q (u,n)) of Y(u,n) and their adjoints since 

T P:q (u,n) = 2- 1 / 2 (U p , q (u,n) + iU* q (u,n)) 

for any p, q, n, u. They are built from variables Xij(u, n) and Xij(u, n), where 

ReX ij:j (u, n) = 2~ 1/2 (ReY i:j (u, n) + ReY j:i (u, n)) 
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lmX id (u, n) = 2~ 1/2 (lmY id (u, n) - ImYj^u, n)) 

ReXj j(u, n) = 2~ 1//2 (ImFj j(u, n) + lmYj :i (u, n)) 

lmXij(u, n) = 2~ 1//2 (ReY^tw, n) — ReYij(u, n)) , 

for any (i,j) e N Ptq and any u s U. These satisfy the assumptions of Theorem 5.1. 
In particular, they are independent due to the fact that for fixed i j and u, n, 
the pairs {ReY it j(u, n), ReYj^u, n)} and {Iml^j(w, n), ImYj-^s, n)} are identically dis- 
tributed since the variance matrices (v Ptq (u)) are symmetric by assumption. Denote 

U Piq (u,n) = Z Ptq (2u- l,n) and U p , q (u, n) = Z Piq (2u, n) 

for any p,q s [r] and us U. Using Theorem 5.1, we can express the asymptotic mixed 
moments in the blocks (Z P)q (u, n)) in terms of mixed moments in symmetrized Gaussian 
operators (u) Ptq (u)), where u s [2t], namely 

^r q (n)(Z puqi ( Ul ,n) . . . Z Pmjqm (u m ,n)) = V q (uj puqi ( Ul ) . . . u} Pm , qm {u m )) 

for any (pi, q±), . . . , (p m , q m ) e J, q e [r], and ui, . . .u m e [2t\. Each of these arrays has 
semicircle distributions on the diagonal and Bernoulli distributions elsewhere. These 
arrays are independent in the sense discussed in Section 3. Now, the linear map 

X : M 2 *(C) - M 2 '(C) 

such that 

X (e p , q {2u-1)) = 2- l / 2 (e Piq (2u-l) + ie Piq (2u)) 
X (e p , q (2u)) = 2- 1 /\e p>q (2u-l)-ie Ptq (2u)) 
for any p,q s [r] and u s [t], induces a unitary map 

Mix) ■ M{U)^M{U) 

where M {%) is the matricially free Fock space of tracial type over the array H = ('Hp,q) 
of Hilbert spaces 

2t 

K P , q = ©C M (tt), 

u=l 

where p, q e [r] , for which it holds that 

M( x )v q = n q , 

M(x)p P , q (^ l)M(x)* = 2- l ' 2 {pU2u 1) + ip p , q (2u)), 
M(x)Pp, q (2u)M( X y = 2- 1 / 2 (p M (2u-l)- i p p>q (2u)) 
for any p,q s [r] and me [t] . Consequently, 

M( X )VpMM(xT = 2- 1/2 (cu Pl9 (2 M - 1) + iu p , q (2u)) 

for any p,q s [r] and us [t]. Therefore, the mixed *-moments in the operators 

Cj m {2u - 1) + ico Pjq (2u) 

can be expressed as the corresponding mixed *-moments in the operators rj Piq (u), where 
p, q s [r] and us [t] , respectively. This completes the proof. ■ 
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Remark 9.1. Let us observe that the family of arrays 

{[C(u)] : u 6 u], 
where each entry of [C( M )] is an operator of the form 

C P ,q( u ) = Z> P ,q( 2u ~ !) + ^p, 9 (2m), 
for any (p, q) e J , is a matrix- valued analog of the family of Voiculescu's circular 
operators [32] related to the circular law. The origins of the circular law in random 
matrix theory date back to the work of Ginibre [14], whereas the papers of Bai [3] and 
Tao and Vu [31] showed its universality (see also the survey paper of Bordenave and 
Chafai [9]). In particular, it is easy to see that each operator ( Ptq (u) is of the form a + ib, 
where a and b are free with respect to the corresponding state ty q and have semicircle 
distributions under ^ q . If p = q, the distributions of o and b are identical, whereas 
if p ^ q, they are identical if and only if d p = d q . The arrays [C(w)] play the role of 
matricial circular operators. 

Theorem 9.1 can be applied to the study of asymptotic distributions of Wishart 
matrices and related products of complex rectangular Gaussian random matrices. For 
that purpose, we shall give a combinatorial formula for the mixed *-moments of the 
operators from the arrays (i] Pjq (u)), where u e 11 (in combinatorics, these operators 
are more convenient than those from {Cp,q{ u )))- First, however, let us present a simple 
example with explicit computations of such *-moments. 

Example 9.1. For simplicity, assume that t = 1 and denote r) Piq = rj Piq (l) for any 
p, q 6 [r]. Of course, 

Vp,q = PpA 1 ) + P*p,q( 2 ) aIld Vp,q = Pp,q(^) + PpA 2 ) 

For any q s [r], we have 

= **(Pm(2)|W2)) = W2) 

yq((Vp,qV;,q) 2 ) = *q(p;, q (2)p p , q (2)p;, q (2) Pptq (2)) 
+ ^q(pl, q (2)pl p {l)p q , P {l)p p , q {2)) 

= b 2 Ptq (2) + W2) Wl) 
^((W;,/) = *,(^(2)^(2)^(2)^(2)^(2)p M (2)) 
+ ^(^(2)^(2)^(2)^ p (l)^ )P (l)p Pi9 (2)) 

+ * ff (p; ff (2) P ; iP (i)p w (i) Pi , ff (2)p; ff (2) ftMI (2)) 

+ ^(^(2)p*, p (l)^,p(l)p* )P (l)^,p(l)p Pl9 (2)) 
+ ^(^(2) P ; p (l)^(2)^(2)^ p (l)p Pi9 (2)) 

= 6^(2) + 3^(2)6 W (1) + W2)&L(1) 

and thus the summands correspond to non-crossing colored partitions. For instance, 
the summands in the last moment correspond to suitably colored partitions shown in 
Example 5.4. In the random matrix setting of Theorem 9.1, if all block variances are 
set to one, the expressions on the right-hand sides reduce to polynomials in d p , d q of 
the form 

d'P'j dp dpdq j *^dp ~~\~ d-pdq 

respectively, since b P)q {k) = d p and b qtP {k) = d q for k e {1,2} (cf., for instance, [16]). 
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The above example shows the connection between mixed *-moments of the considered 
operators and non-crossing colored partitions. The main feature of the combinatorics 
for the non-Hermitian case is that in order to get a pairing between two operators, one 
of them has to be starred and one unstarred. This leads to the following definition and, 
consequently, to an analog of Proposition 4.3. 

Definition 9.1. We say that n e MC 2 m is adapted to ((zt>i, u±, ei), . . . , («/ m , u m , e m )), 
where zv k = {p k ,q k } and (p k ,q k ,u k ) e [r] x [r] x [t] and e k e {1,*} for any k, if 
there exists a tuple ((^, iti), . . . , (v m , u m )), where v k e {(p k , q k ), (q k ,p k )} for any k, 
to which ir is adapted and €j ^ e k whenever {j, k} is a block. The set of such 
partitions will be denoted by jVC^((o^, ui, ei), . . . , («/ m , u m , e m )). Its subset consist- 
ing of those partitions which are in jVC^ 9 ((zyi, u±), . . . , {rv m , u m )) will be denoted 

Proposition 9.1. For any tup/e Mi, ei), . . . , (zv m , u m , e m )) and q e [r], m e N, 
where iv k = {p k , q k } and p k , q k e [r], u k e [t] and e k e {1, *} for each k, it holds that 

* 9 (C,^)---C,*>-)) = S ^(^z) 

7rejVC^ lg ((7yi,ui,ei),...,(w m ,u m ,e m )) 

where f is the coloring of n defined by ((zi\,Ui), . . . , (zv m ,u m )) and by the number q. 

Proof. The proof is similar to that of Proposition 4.3. As in the case of u) Ptq (u), 
each operator rj Piq (u) and rj* Ju) is a sum of a symmetrized creation operator and a 
symmetrized annihilation operator. However, in the non-self-adjoint case, these sum- 
mands are labelled by two different numbers, 2u— 1 and 2w, respectively. Therefore, the 
non-crossing pairings can be of two types: (r) Ptq (u), r)* q (u)) or (r)* q (u), r) Pjq (u)), which 
means that the associated non-crossing partitions must satisfy the additional condition 
€j e k whenever {j, k} is a block. This completes the proof. ■ 



10. WlSHART MATRICES AND RELATED PRODUCTS 

The main point of this section is to show that the framework of matricial freeness is 
general enough to include sums and products of independent Gaussian random matrices 
with block-identical variances in one unified scheme, although the formalism is much 
more general and should give other interesting examples. 

In particular, we can reproduce certain results concerning products B{n) of rectan- 
gular Gaussian random matrices and the asymptotic distributions of B{n)B*{n) under 
the trace composed with classical expectation as n — > go. In particular, if B{n) is just 
one Gaussian random matrix, the matrix 

W{n) = B(n)B*(n) 

is the complex Wishart matrix [37]. The limit distribution of a sequence of such matrices 
is the Marchenko-Pastur distribution (which also plays the role of the free Poisson 
distribution) with moments given by Catalan numbers. 

The original result on the asymptotic distributions of the Wishart matrices is due to 
Marchenko and Pastur [25], but many different proofs have been given (see, for instance 
[15,27,35]). In the case when B{n) is a power of a square random matrix, it has recently 
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been shown by Alexeev et al [1] that the limit moments are given by Fuss-Catalan num- 
bers (products of independent random matrices have also been studied recently [2,10]). 
The distributions defined by the Fuss-Catalan numbers were explicitly determined by 
Penson and Zyczkowski [29]. A random matrix model based on a family of distributions 
called free Bessel laws constructed from the Marchenko-Pastur distribution by means 
of free convolutions was given by Banica et al [4]. In turn, asymptotic freeness of inde- 
pendent Wishart matrices was proved by Voiculescu [32] (see also the stronger versions 
of Capitaine and Casalis [11] and Hiai and Petz [18]). Let us also mention that explicit 
formulas for moments of Wishart matrices were given, for instance, by Hanlon et al [17] 
and by Haagerup and Thorbjornsen [16]. 

In our framework, it suffices to use one symmetric off-diagonal (rectangular or square, 
balanced or unbalanced) block to study the asymptotic distribution of the Wishart 
matrix. More generally, using the same operator algebras as in Section 5, we can also 
study the asymptotic distributions of the matrices W(n) in the case when each B(n) is 
a sum of independent rectangular Gaussian random matrices as in [6], namely 

B(n) = Y(u u n) + Y(u 2 , n) + . . . + Y(u m , n), 

assuming they have the same dimensions, as well as in the case when B{n) is a product 
of independent rectangular Gaussian random matrices 

B(n) = Y(u 1 ,n)Y(u 2 , n) . . . Y(u k , n), 

assuming their dimensions are such that the products are well-defined. Using Theorems 
5.1 and 9.1, one can study non-Hermitian Wishart matrices as well. In both cases, it 
suffices to take a suitable sequence of the off-diagonal symmetric blocks of a family 
of independent Gaussian random matrices (the first case) or even one Gaussian ran- 
dom matrix matrix (the second case). Since we know how to compute the asymptotic 
joint distributions of symmetric blocks, we can immediately obtain results concerning 
asymptotic distributions of such products. 

In our treatment of products, it will be convenient to use sets AfC^W^p) of non- 
crossing W^p-pairings, where W kjP is a word of the form 

W k>p = (12... pp* ...2*l*) k 

where m = 2kp, by which we understand the set of non-crossing pair partitions of the 
set [m] in which all blocks are associated with pairs of letters of the form {j,j*}. If 
7r e AfC 2 m (W k , p ), we denote by TZ^tx) and TZ*(tt) the sets of the right legs of it which 
are associated with j and j*, respectively. 

Definition 10.1. Define homogenous polynomials in variables d±, d 2 , ■ ■ ■ , d p+ i of the 
form 

Pk(di, d 2 , ■ ■ ■ , d p+ i) = ^ dl^dl, 2 ^ . . . (ip+x 1 ^ 

neUCl kp {W Kp ) 

for any k,p e N, where 

r,(7r)H^MI + l^*-iMI 
for any tt g MC\ kp iyV) and 1 ^ j ^ p + 1, where we set TZq(tt) = and 7^ p+ i(7r) = 0. 

Our random matrix model for products of independent Gaussian random matrices 
will be based on products of a chain of symmetric random blocks 

Tjj +1 (n) = Tjj +1 (u, n), where j e [p] and n g N, 
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with fixed us U and thus omitted in our notations (taking different it's does not change 
the computations since these blocks ore independent anyway). These blocks are em- 
bedded in a Gaussian random matrix Y{n) with asymptotic dimensions d\, d 2 , ■ ■ ■ , d p+ \ 
which will be the variables of our polynomials, namely 

Y{n) = j]T h3+1 {n) 

where each Tjj+i(n) consists of two blocks, Sjj + i(n) and Sj + ij(n), placed symmet- 
rically with respect to the main diagonal. The blocks 7} J+ i(n) are not assumed to 
be Hermitian, thus we shall use Theorem 9.1 when discussing the limit distribution 
of W(n). Using products of Hermitian random blocks is also possible, especially in 
the simplest case of Wishart matrices. However, it is more convenient to use products 
of non-Hermitian symmetric random blocks to give a combinatorial description of the 
limit moments of W(n) in the case when B(n) is a product of order higher than two. 

In the proof given below, we use ideas on the enumeration of non-crossing pairings 
given by Kemp and Speicher [20] (related to the results in the book of Nica and Speicher 
[26]). 

Theorem 10.1. Under the assumptions of Theorem 9.1, suppose that 11 consists of 
one element and Vjj+i = 1 for j e [p], where p is a fixed natural number, and let 

B(n) = T^ 2 {n)T 2 ,z{n) . . .T P:P+1 (n) 

for any n s N. Then, for any k e N, 

lim n(n) ((B(n)B*(n)) k ) = P k (d 1} d 2 , . . . , d p+1 ) 
where d 1 ,d 2 , ■ ■ ■ , d p+ i are asymptotic dimensions. 
Proof. By Theorem 9.1, we obtain 

lim T l {n){{B{n)B*{n)f) = ^i((^*) fe ) 

n^oo 

for any natural k, where 

V = Vl,2V2,3 ■ ■ ■ Vp,p+i 

for any given p. In view of Proposition 9.1, we can express the right-hand side in terms 
of non-crossing partitions 

TTeNC 2 ml ((wi,t±),...,{u> m ,£ m )) 

where m = 2kp and 

{iv 1+2jp , . . . , w p+2jp ) = ({l,2},...,{p,p+ 1}) 

(^ P +i+2jp, • • • , w 2p+2jp ) = ({p,p+ 1}, . . . , {1,2}) 

for any < j < k — 1. In the notation for non-crossing pair partitions adapted to 
((zi>i,ei), . . . ,(zi/ m ,e m )) we omit u\,...,u m since in the considered case we have one 
array (rfij). Each is a sum of the form 

Vjj+i = Pjj+iiX) + P*j+i( 2 )' 
but in the computation of moments, the information about which creation-annihilation 
pairs (labelled by 1 or 2) contribute to the moment ^i((r]r]*) k ) is encoded in e = 
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(ei, . . . , e m ). Namely, the pairs (p* J+1 (l), p jJ+1 (l)) or (p* J+1 (2), p jJ+1 (2)) contribute a 
number from the set {fyj+i, kj+ij} = {dj, dj+i} depending on whether (Vjj+n Vj,j+i) or 
(r]jj + i, Tjjj+i) is associated with the given block of it. Moreover, only operators with the 
same matricial indices can form a pairing and thus HC 2 m i((a>L, ei), • • • , (aw, e m )) can be 
put in one-to-one correspondence with the set NC\ kp {Wk,p)- If now it e NC\ kp {Wk,p)-, 
then we will identify bi(n, /) with the corresponding expression of Proposition 9.1. 
Recall that 

kp 

b l (nJ) = Y\b l (n 3 J) 

3 = 1 

where {tt\, . . . , 7Tk P } are blocks of it and where the coloring / is uniquely determined by 
7r and by the coloring of the imaginary block (in this case, the latter is equal to one 
by our choice of T\(n)). Since the block variances are set to one, we can express each 
bi(n, f) as a monomial of order m = kp in asymptotic dimensions of the form 

where the exponents are natural numbers which depend on it and for which 

ri(yr) + r 2 (7r) + . . . + r p+ i(7r) = kp. 

Let us find an explicit combinatorial formula for these numbers. When computing 
bi(7Tk, f) for various blocks of it, we can see that asymptotic dimensions are assigned to 
the right legs associated with letters {j,j* : j e [p]} according to the rules 

j — > dj and j* — > 

for any j s [p] . This leads to the formula 

rj (n) = |^(tt)| + I^MI for j e [p + 1], 

where it e NC 2 m (W), and thus 

*i«w*) fc )= 2 < 2W < 2W ...<ri lW , 

which completes the proof. ■ 

The special case of p = 1 corresponds to Wishart matrices and the Marchenko-Pastur 
distribution [25] with shape parameter t > 0, namely 

a/ (x — a) (b — x) , . , 
ft = max{l - t, 0}S + V V 2 J^ J -l [aM {x)dx 

where a = (1 — Vt) 2 and b = (1 + V^) 2 (see Corollary 10.1). At the same time, g t is the 
free Poisson law in free probability. It has been shown by Oravecz and Petz [27] that 
its moments are the Narayana polynomials 

k 

N k (t) = Y l N(k,j)t j 

3=1 

for any fceN, with coefficients 

1 (k - 1\ / fc 



■V(/.-..y) , . , M 
3 \3 ~ 1/ VJ 
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called Narayana numbers. These numbers are obtained in several different enumera- 
tions. For instance, N(k,j) is equal to the number of Catalan paths of lenght 2k with 
j peaks. However, we will use another enumeration related to right legs of blocks of the 
associated non-crossing partitions. 

Proposition 10.1. The Narayana number N(k,j) is equal to the number of those non- 
crossing pair partitions ty e MC\ k which have j even numbers in the set 1Z(n). 

Proof. Our proof is based on the enumeration derived by Osborn [28], in which 
N(k,j) is the number of Catalan paths of lenght 2k with 2j edges lying in odd bands, 
where an odd band is a set of the form Ix [i, i + 1] in the plane M 2 , where i is even, in 
the standard setting, in which a Catalan path begins at (0, 0) and ends at (2k, 0). Now, 
let us associate with each Catalan path a non-crossing pair partition -k g NC\ k in the 
canonical way. Then, if {/, r} is a block of n with / < r, where r is even, then the edges 
corresponding to I and r lie in the odd band. Similarly, if r is odd, then these edges 
lie in the even band. Consequently, N(k,j) is equal to the number of non-crossing pair 
partitions of [2k] having j right legs which are even and k — j right legs which are odd. 
This completes the proof. ■ 

Example 10.1. In the path shown below all edges labelled by a lie in the odd bands 
and all edges labelled by b lie in the even bands. Therefore, there are 6 edges lying in 
the odd bands and 8 edges lying in the even bands, thus j = 3. 




Corollary 10.1. If p = 1 and di > 0, then under the assumptions of Theorem 10.1 we 
obtain 

lim n(n) ((B(n)B*(n)) k ) = d^N^/d,), 

n—>cc \ / 

and thus the limit distribution is the di- dilation of the Marchenko-Pastur distribution 
Qt with the shape parameter t = d<ijd\, denoted Qd 2 ,d x - 

Proof. The letters 1 and 1* correspond to odd and even numbers, respectively, if we 
use the notation as in the proof of Theorem 10.1. In turn, the number of non-crossing 
pair partitions which have j right legs labelled by 1* is given by N(k,j) by Proposition 
9.1 and thus 

k 

Pk(dud 2 ) = '£ i N(k,j)4<%- j = dtN^/dJ 
j'=i 

which proves the first assertion. The second one is clear since the transformation on 
moments m k — > \ k m k leads to the dilation of the corresponding measures. ■ 



Remark 10.1. The usual settings for Wishart matrices are often slightly different since 
different normalizations are used [10,19]. For instance, if B(N) e GRM(m(N), N, l/N) 
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is a Gaussian random matrix of dimension m(N) x N and 1/N is the variance of each 
entry, where N e N, it is assumed that 

m(N) 
t = lim ——^ 

JV->oo iV 

and then one computes the limit distribution of B*(N)B(N) under normalized trace tr^r 
composed with classical expectation. One obtains Narayana polynomials in t as the limit 
moments, with the corresponding Marchenko-Pastur distribution with shape parameter 
t. In our model, we embed B{n) and B*{n) in a larger square matrix of dimension n, 
in which B(n) = S\ :2 (n) and B*(n) = S* 2 (n) = S 2: i(n) are off-diagonal blocks and we 
compute the limit moments of B(n)B*(n) under Ti(n). In order to directly compare 
these two approaches, we set iV = ni{n) and m(N) = n 2 {n) to get t = d 2 /d\ and, since 
the variance in our approach is 1/n and in Marchenko-Pastur's theorem the variance 
is 1/N, the kth moment in the limit moment obtained in Marchenko-Pastur's theorem 
must be multiplied by d\ to give our asymptotic product polynomial. 

Example 10.2. If p = 1, the lowest order polynomials of Definition 10.1 are 

Pi(d 1 ,d 2 ) = d 2 

P 2 (di,d 2 ) = d 2 di + d\ 

P 3 (d 1 , d 2 ) = d 2 d\ + t>d\d x + d\ 

P 4 (d u d 2 ) = d 2 d\ + Qd\d\ + Qd\d x + d\ 

and can be obtained directly from Corollary 10.1 since the corresponding Narayana 
polynomials are N x {t) = t, N 2 (t) =t + t 2 , N 3 (t) =t + 3t 2 + 1 3 , N 4 (t) =t + 6t 2 + 6t 3 + 1 4 , 
respectively. 

Corollary 10.2. If d± = d 2 = . . . = d p+ i = d, then 

P k (d 1 ,d 2 ,...,d p+1 ) = d k PF(p,k) 

where F(p,k) = ^p[(^ k ^ k ) are Fuss- Catalan numbers. 

Proof. If di = ... = d p+ i = d, then Pk(d\,d 2 , . . . ,d p+ \) is equal to d kp multiplied 
by the number of non-crossing pair partitions adapted to the word Wk, P by Theorem 
10.1. It is well-known [20] that the latter is equal to the Fuss-Catalan number F(p, k). ■ 

Example 10.3. If p = 2, we obtain by Theorem 10.1 the limit moments 

^i((w*) fc ) = p k(di, d 2 , d 3 ), where r] = rj 1:2 rj 2:3 . 
For instance, contributions to P 2 (d±, d 2 , d 3 ) are 

^l(p2,lP3,2p3,2p2,lp2,lP3,2p3,2p2,l) = d 2 d 3 
^l(p2,lP3,2p3,2p* 2 Pl,2p3,2p3,2p2,l) = ^1^2^ 
^l(p2,lP3,2p2,3Pl,2pl,2p2,3p3,2p2,l) = ^1^3 

with the corresponding partitions of the associated word W 2y2 = 122*1*122*1* 
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I n I I n I Innnl ml 

respectively. Thus, 

P 2 (di, d 2 , d 3 ) = d 2 d^ + d\d 2 dl + did%d 3 . 

If we define multivariate generalizations of Narayana polynomials N k (ti,t 2 , . . . ,t p ) by 
the formula 

P k (di, d 2 ,..., dp+i) = d k 1 p N k (d 2 /d 1 , d 3 /d 1 , dp+x/dx), 

we obtain 

N 2 (h,t 2 ) = t\t\ + tit 2 2 + t\t 2 . 

These polynomials and their coefficients are of some interest from the combinatorial 
point of view and will be studied elsewhere in more detail. 

The asymptotics of the Wishart matrix can be generalized in another direction, 
namely we can study the asymptotic distribution of W(n) = B(n)B*(n), where B(n) 
is now a sum of independent rectangular Gaussian random matrices. In our theory, we 
can compute the corresponding limit moments under the trace by choosing each of these 
matrices to be the rectangular off-diagonal blocks Si t2 (u, n) and using either Hermitian 
or non-Hermitian symmetric blocks containing them. In the proof of the theorem given 
below, we use Hermitian symmetric blocks. 

Theorem 10.2. If {S(u,n) : u e U) is a family of rectangular independent complex 
Gaussian random matrices with identically distributed entries for any ?ieN and 

B{n) = Y J S{u,n), 

ueu 

then the limit distribution ofW{n) = B{n)B*{n) under the partial trace T\{n) is Qd 2 ,d 1 - 

Proof. Let us treat each matrix S(u, n) as the block Si t2 (u, n) in a Hermitian random 
matrix Y(u,n). First, let us show that 

lim r 1 (n)(( J B(n) J B*(n)) fe ) = *i((£ li2 ) 2fc ) 

n^co 

for any k e N for some Ui t2 = u>i t2 + ui 2t \. If we embed the blocks Si t2 (u,n) is the 
Hermitian blocks T li2 (u,n), then 

ri(n)(S lt2 (u 1 ,n)Sl 2 (u 2 ,n) . . . S lt2 (u 2k - 1 ,n)Sl 2 (u 2k ,n)) 

= Ti(n)(T lj2 (ui,n)Ti t2 (u2,n) . . . T lj2 (u 2 k-i, n)T h2 (u 2k , n)) 

for any u±, . . . u 2k e U and n 6 N since each Hermitian block Ti t2 (uj, n) consists of two 
blocks, Si i2 (uj,n) and its adjoint S* 2 (uj,n), and the partial trace Ti(n) selects only 
vectors from the set {e^ : i e N\} onto which the product acts. Therefore, by Theorem 
5.1, we obtain 

lim r 1 (n)((B(n)B*(n)) k ) = *i((£i, 2 ) 2fc ) 

where 

ueu 
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and cDi )2 = Wi j2 + 0*2,1, with o>i, 2 and 0*2,1 associated with numbers 

h,2 = 2j & i,2(m) and 6 2 ,i = 2] fe 2,i(«)- 

In the computation of moments of 0*1,2 under as in Example 5.2, these numbers are 
assigned to blocks of even and odd depths, respectively. Similar arguments to those in 
the proofs of Proposition 10.1 and Corollary 10.1 allow us to conclude that the moments 
of u\ 2 under \E*i agree with the moments of the e?i-dilation of the Marchenko-Pastur 
distribution with the shape parameter t = 62,1/61,2 = ^2/^1 (by the definition of Her- 
mitian blocks, we must have 1*1,2 (u) = 1*2,1 (u) for any ue 11). ■ 

Example 10.4. The above limit distributions are closely related to those given by the 
Double Convolution Formula. If J7 = {(1, 2), (2, 1)}, then the symmetrization of the 
limit distribution of W(n) under T\(n) is the distribution of the form 

A* = \*ljejt*j = «2,1 ffl «1,2 

corresponding to \E*i, where fij is the Bernoulli law concentrated at iy/bj for j e J ' . 
By the symmetrization of a probability measure v on [0, 00) with moments v{n) we 
understand the symmetric probability measure z/ sym on R with even moments z/ sym (2n) = 
z/(n) and odd moments equal to zero. An analogous result is obtained for the limit 
distribution under r 2 (n) corresponding to \E* 2 - 

Remark 10.2. Alternatively, we can treat each matrix S(u,n) as the block S 1j2 (m, n) 
in a non-Hermitian random matrix Y(u, n). Then again we can take 

B(n) = ^T li2 («, n ), 

ueu 

except that symmetric blocks Ti, 2 (it,n) are not Hermitian. Using Theorem 9.1, we 
obtain 

hm r 1 (n)((B(n)S*(n)) fc ) = 9^^)"), 

n^co ' 

for any k e N u {0}, where 

ueu 

which gives the same limit distribution as in Theorem 10.2. 

Theorem 10.3. Under the assumptions of Theorem 10.1, the family of Wishart ma- 
trices 

W(n) = {S(u, n)S*(u, n) : u e U] 
is aymptotically free under Ti(n). 

Proof. As in the proof of Theorem 10.2, the Wishart matrices will be built from 
blocks 5*i i2 (m) and their adjoints, 

W(u,n) = Si i 2(u,n)S* 2 ( u , n ) 

for any u e U and n e N. Using Theorem 5.1 and arguments similar to those in the 
proof of Theorem 10.2, we can replace Si i2 (-u,n) by Hermitian T 12 (-u,n), which gives 

lim Ti(n)(W(ui,n) . ..W(u k ,n)) = 9 1 {ul 2 {u 1 ) . ..ul 2 (u k )) 
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for any ui, . . .u k e U and n e N. It suffices to prove that the family {uj\ 2 {u) : u e U} 
is free with respect to More generally, we will show that the family of algebras 
{A(u) : us <U), where 

A(u) = C<a)i j2 (u), r-1,2, r 2 ,i> 

and ri,2, r 2 ,i are given by Definition 2.5 (we have li,2 = ri, 2 + r 2 ,i), is free with respect 
to Any polynomial from A(u) is a linear combination of ri,2, r 2 ,i and even powers 
of cDi i2 (w) interchanged with r 12 and r 2 ,i. However, using the relations 

Pl,2( u )pi,2( u ) = b i,2(u)r h2 , pl 1 (u)p 2 ,i(u) = b 2 , 1 (u)r 2A , 

and 

ri,2pi,2(u) = 0, r 2 ,ipi, 2 (w) = pi,2(u), r 2 ,ip2,i(u) = 0, r 1)2 p2,i{u) = p2,iW, 

we can observe that any polynomial a m e A(u m ) n Ker^i, when acting onto Q±, reduces 
to a linear combination of products of the form 

(pl,2(u m )p 2A (u m )) n 

where n > 0, since the remaining terms give zero (in particular, ri, 2 f2i = 0). The 
same is true when we then act succesively with polynomials a m _i, . . . ,a±, where a« e 
A(ui) n Ker^! and m ■ ■ ■ ¥= u rn . We obtain a linear combination of vectors of the 
form 

(ei, 2 (ui) ® e 2:1 ( Ul )) ni ® . . . (x) (e 1)2 (7i m ) ® e 2jl (M m )) nm 

where n\ + . . . + n m > 0, and thus it is orthogonal to Qi, which gives the freeness 
condition 

*i(aia 2 ...a m )=0. 
It remains to observe that li,2 is the projection onto the Hilbert space 

o)i e cn 2 e ® M,fc), 
k 

which is left invariant by all elements of each A(u), which completes the proof. ■ 

Let us finally establish a relation between the limit distributions of Theorem 10.1 
and free Bessel laws 7r Pi t of Banica et al [4] expressed in terms of the free multiplicative 
convolution 

where n = Q\ is the standard Marchenko-Pastur distribution with the shape parameter 
equal to one. Random matrix models for these laws given in [4] were based on the 
multiplication of independent Gaussian random matrices. 

Let n s = U s fi be the probability measure on the real line defined by 

for any s > 0, where G^ s and are Cauchy transforms of fi s and fj,, respectively. 

Theorem 10.4. Ifdi, d 2 , . . . , d p+ \ > andp > 1, then the limit distribution of Theorem 
10.1 takes the form 

f*i = U s (g dud2 g d2jd3 ... Q dp+1 ,d P ) 

where s = d p /d\. 



LIMIT DISTRIBUTIONS OF RANDOM MATRICES 61 

Proof. We have 

Ti(ra) ((Si >2 (rO . . -^+1(71)5^+! (ra) . . . S^fa))*) 

= na/m • r 2 (n) ((S 2 » . . . ^ +1 (n)S* p+1 (n) . . . ^>)S? >2 (n)Si >2 («)) fc ) 

= n 2 /ni • r 2 (n) ((T 2>3 (n) . . . T P:P+1 (n)T P:P+1 (n) . . . T 2 , 3 (n)T 1>2 (n)T 1;2 (n)) k ^ 
which tends to 

d 2 /d 1 -^ 2 {{uo 2 ujl 2 ) k ), 

as n — > 00, where cD = 0)2,3 • • • &p,p+i- Now, the pair {u> 2 , ui\ 2 } is free with respect to ty 2 . 
The proof of this fact is similar to the proof of the asymptotic freeness of independent 
Wishart matrices given in Theorem 10.3 and is omitted. Now, the ^-distribution of 
uj\ 2 is Qd u d 2 by Corollary 10.1. Therefore, 

= U d2 / dl (g dud2 @^ 2 ) 

where fi 2 is the ^-distribution of cD 2 . If ui is a product of at least two operators, the 
same procedure is applied to fi 2 , which gives 

^2 = U ds / d2 (g d2td3 S/x 3 ) 

where /x 3 is the ^-distribution of (0)3,4 . . . cD PiP+ i) 2 . We continue this inductive procedure 
and observe that the last step gives 

fJ>p-i = U dp / dp _ 1 (g dp _ 1)dp \x\ g dp+ltdp ), 

which gives the desired formula for Hi since the operation U s commutes with \x\ and we 
can collect all operations of this type to get U d2 / dl . . . U dp i dp _ x = U dp/dl . ■ 

Corollary 10.3. In particular, if d\ = d 2 = . . . = d p = 1 and d p+ i = t, the limit 
distribution of Theorem 10.1 is the free Bessel law 7i Ptt . 

Proof. This is an immediate consequence of Theorem 10.4. ■ 

11. Appendix 

The definition of the symmetrically matricially free array of units given in [23] should 
be strenghtened in order that the symmetrized Gaussian operators be symmetrically 
matricially free. This requires certain changes to be made in [23] which are listed below. 

(1) Condition (2) of Definition 8.1 in [23] should be replaced by condition (2) of 
Definition 3.4 given in this paper. For that purpose, one needs to distinguish 
even and odd elements. 

(2) The proof of symmetric matricial freeness of the array of symmetrized Gaussian 
operators was stated in Proposition 8.1 in [23] and the proof was (unfortunately) 
omitted. Using Definition 3.4, this fact is proved here in Proposition 3.5. 

(3) Definition 2.3 in [23] should be replaced by Definition 4.3 given in this paper. 
In principle, it is possible to modify the old definition and use conditions on 
intersections of unordered pairs, but it is much more convenient and simpler to 
phrase the new definition using sequences of ordered pairs. 
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(4) The proof of Theorem 9.1 in [23] needs to be slightly strengthened since the 
classes NC? m q {{pi, <?i}, • • • , {p m , Qm}) can be smaller, in general, than those con- 
sidered in [23] since Definition 4.3 is stronger than Definition 2.3 in [23]. There- 
fore, we need to justify that if 71(7) e NC 2 mq \NC 2 m>q {{pi, qi},..., {p m , q m }), then 
the corresponding mixed moment of symmetric blocks under r q (n) vanishes. 
Clearly, in order that this moment be non-zero, the imaginary block of 7r(7) 
must be colored by q. Then all blocks of 71(7) of depth zero must be labelled by 
pairs {pi,qi} which contain q in order that the corresponding symmetric blocks 
T p . >q .(n) act non-trivially on vectors ej,j s N q . Supposing that qi = q for all 
such pairs {pi,qi} assigned to the right legs of n^j), we set r% = (f>i, <??,)• Then, 
in turn, all symmetric blocks which correspond to blocks of 71(7) of depth one 
whose nearest outer block is labelled by given {pi, q^} must be labelled by pairs 
{pk,qk} which contain pi. Supposing that g& = pi for all such pairs {pu,qk} as- 
signed to the right legs of 71(7), we set Vk = (pk,qk)- We continue in this fashion 
until all symmetric blocks and the corresponding blocks of 71(7) are taken into 
account. Finally, if {z, j} is a block, where i < j and Vj = (pj, qj), then we must 
have Vi = (qj,Pj) in order that the action of T PiM (n) be non-trivial, which follows 
from an inductive argument starting from the deepest blocks. Consequently, in 
order to get a non-zero contribution from the corresponding mixed moment of 
symmetric blocks, which is proportional to the product of variances as in the 
proof in [23], there must exist a tuple (vi, . . . , v m ) such that Vi s {(p^, qi), (qi,Pi)} 
to which 71(7) is adapted. Therefore, the conditions of Definition 4.3 are satis- 
fied. 
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